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PREFACE 


The  Twelfth  Annual  Space  Surveillance  Workshop  sponsored  by  ESC  and  hosted  by  MIT 
Lincoln  Laboratory  was  held  on  5,  6  and  7  April  1994.  The  purpose  of  this  series  of  workshops  is 
to  provide  a  forum  for  the  presentation  and  discussion  of  space  surveillance  issues. 

This  Proceedings  (Volume  II)  contains  unclassified  presentations  fi'om  this  workshop 
received  after  the  deadline  for  Volume  I.  The  papers  were  reproduced  directly  from  copies 
supplied  by  their  authors  (with  minor  mechanical  changes  where  necessary).  It  is  hoped  that  this 
publication  will  enhance  the  utility  of  the  workshop. 


Kurt  P.  Schwan 
Editor 


The  Quest  for  Russian  Optical 
Surveillance  Technolosv 
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Technology  of  interest  -  -  Optical  Phase  Conjugation  (OPC) 


Illustration  of  the  Aberration-Correction  Properties 

of  Phase  Conjugation 


Comparison  of  a  Conventional  Mirror  Reflection 
and  a  Phase-Conjugate  Mirror  Reflection 
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Optical  Phase  Conjugation  by  Four-Wave  Mixing 


Potential  Applications  of  OPC 
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Space  debris  tracking 


Space  Catalog  Snapshot  (22Feb) 


This  does  nfit  reflect  the  total  number  of  earth-orbiting  objects! 


Application  of  Interest: 
Snacetrack  of  Small  Obiects 
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1  km  positional  accuracy  (1  sigma) 

0.25  m/sec  orbital  speed  accuracy  ( 1  sigma) 


Russian  Institutes  Visited 
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fassat  ^JMznny-JNovgoroa; 
Laser  Systems,  Inc.  (Moscow) 


Team  Observations  (trip  1) 
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We  sensed  a  real  willingness  to  show  us  more  in  time 


Team  Observations  (trip  2) 
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Promising  areas  in  satellite  catalog  maintenance 

Invitations  to  visit  optical  facilities  in  Caucasus  and  Central  Asia 


A  Quick  Tour  of  Russia 


C/D 

cd  C 

^  s 

^  o 

S  CO 

°  S 

<T3 
.  C 

cd  *cr3 

•1-H 

CO  H 
CO  ^ 

o  S 

W)§ 

CX 

2  ’S  <D 

§  S  ^ 

H  S 

Sj3 


a  - 


3  T3 

o  c 

>.  13 

<U  4^ 

> 

<U 

a 
<u 
413 


>> 

.  § 
0) 

s  > 

s 

S  o 

V->  4-> 


CD 

c 


3  On 

cl  ^ 

W  00 

^  T-H 

^  e 

CO 

*3 

ar^‘3 

v-i  CO 

cd  3 

Soi 


General  Observations 
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General  Observations 


General  Observations 


Utility  of  the  MSX  UVISI  Instrument  to  Space  Surveillance: 

Imaging,  Tracking  and  Identification 

P.  F.  Bythrow,  R.  L.  Waddell,  Pat  Murphy  and  J.  F.  Carbary 

The  Johi\s  Hopkins  University  Applied  Physics  Laboratory,  Laurel,  MD  20723 


ABSISAa 

The  UVISI  Instrument  of  the  Midcourse  Space  Experiment  (MSX)  consists  of  four  imagers  and 
five  spectrographic  imagers  (SPIMs)  using  intensified  CCD  detectors  for  imaging  and  determination  of 
target  spectral  content.  UVISI  incorporates  1750A  based  Image  and  Track  processors  to  conduct  target 
selection,  track  file  computations  and  trajectory  determination.  The  UVISI  Imagers  use  a  filter  wheel 
arrangement  to  cover  the  wavelength  domain  from  110  nm  to  900  nm  in  selectable  bands.  Imager 
sensitivity  is  of  the  order  of  5  photons /cm^sec.  Intrascene  dynamic  range  is  4096,  while  an  interscene 
dynamic  range  of  -10^  is  allowed  by  the  micro  chaniwl  plate  (MCP).  This  broad  d)mamic  range  allows 
UVISI  to  view  dim  objects  in  deep  space  as  well  as  near  earth  objects  agamst  a  sunlit  earth  background. 
SPIMs  cover  the  same  wavelengto  regime  as  the  imagers  and  have  similar  sensitivity.  The  SPIMs  have 
spectral  resolution  of  -0.5  -  4.0  run.  This  paper  will  describe  dtese  characteristics  in  greater  detail  and 
discuss  how  UVISFs  unique  capabilities  can  applied  to  the  space  surveillaitce  task.  Speafically  UVISI 
has  the  potential  for  discrintination,  viewing  and  tracking  Resident  Space  Objects  (RSOs)  against  a  sunlit 
earth  and  earth  limb  as  well  as  deep  space.  UVISI  also  has  the  uitique  potential  to  conduct  multispectral 
Space  Object  Identification  (SOI)  by  the  identification  of  specific  reflectance  spectra  from  various  RSO 
component  materials. 


l.IFTraQDUaiQN 

The  Midcourse  Space  Experiment  (MSX)  is  a  spacecraft  mission  spoitsored  by  the  Ballistic  Missile 
Defense  Organization.  The  mission's  principal  design  objective  is  to  gather  data  on  targets  and 
backgrounds,  both  deep-space  and  near-Earth.  The  spacecraft  has  been  developed  aiKl  integrated  by  the 
Johns  Hopkins  University  Applied  Physics  Laboratory  (APL),  and  is  due  to  be  launched  in  November 
1994.  Major  sensor  systems,  developed  by  APL,  Utah  State  University’s  Space  Dynamics  Labota  *  h  v  and 
MIT s  liiKoin  Laboratory,  include  the  UVISI  ultraviolet  and  visible  imagers  and  spectrographs  (APL),  the 
SPIRIT  III  infrared  imager  and  interferometer  (USU/SDL),  and  the  SBV  space-based  visible  sensor 
(MIT/LL).  The  spacecraft  is  designed  for  a  five  year  h'etime  after  being  launched  into  a  near-sun- 
synchronous  orbit  of  -  9(X)  kilometers,  inclined  about  992  degrees.  A  sophisticated  attitude  control  and 
determination  system  and  a  closed-loop  tracking  system  are  integral  to  foe  data  collection  experiments. 

The  UVISI  instrument  derives  from  a  succession  of  ultraviolet  and  visible  instruments  flown  on 
ofoer  .\PL  orbital  missions  over  foe  past  decade  for  the  Department  of  Defense  aixl  NASA.  The  AIM 
(Auroral  Ionospheric  Mapper)  flew  on  the  HILAT  spacecraft  in  1983^,  and  its  cousin  the  AIRS  (Auroral 
Ionospheric  Remote  Sensor)  flew  on  foe  Polar  BEAR  spacecraft  in  1986^.  Both  instruntents  made 
measurements  of  foe  ultraviolet  aurora  against  a  suirlit  Earth.  The  Hopkins  Ultraviolet  Telescope  made 
astronomical  measurements  of  77  celestial  objects  from  the  Shuttle  STS-35  mission  in  19903-  APL 
instruments  also  flew  on  foe  Department  of  Defense/Strategic  Defense  Initiative  series  of  Delta  missions. 
These  instruments  made  foe  first  exoatmospheric  observations  of  UV  plume  phwomena  and  foe  first 
multi-spectral  observations  of  an  RSO  from  orbit^'^A 
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UVISI  consists  of  a  suite  of  five  spectrographic  imagers,  four  imagers,  an  image  processor 
subsystem,  and  various  other  electronic  control  and  data  handling  subsystems.  Figure  1  shows  the  block 
diagram  of  the  UVISI  system  and  indicates  the  major  subsystems  of  the  instrument.  Each  sensor  (SPIM  or 
Imager)  has  an  independent  power  supply  (converter)  and  sensor  electronics  units  (SEU)  that  allows  it  to 
operate  autonomously  should  any  of  the  other  sensors  fail.  (The  SPIMs  are  numbered  1  through  5 
according  to  the  spectral  regime,  with  SPIMl  sampling  the  far  ultraviolet  and  SPIMS  sampling  the  visible. 
The  imagers  are  labeled  according  to  spectral  regime  and  field  of  view.  Thus,  lUN  indicates  imager, 
ultraviolet,  narrow-field.)  The  data  control  s)^tem  (DCS)  accepts  commands  from  the  spacecraft  and 
formats  and  sends  digital  data  from  die  sensors  back  to  the  spacecraft  and  to  the  UVISI  image  processor. 
The  DCS  has  a  backup  in  case  of  failure  of  the  primary  subsystem.  The  image  processor  reduces  the  raw 
image  data  from  any  one  of  the  four  imagers  and  directs  the  results  to  the  spacecraft  tracking  processor. 
The  track  processor  may  then  use  target  data  from  the  image  processor  to  slew  the  spacecraft  or  the 
spacecraft  may  be  pointed  open  loop  and  the  track  file  accumulated  and  stored  for  later  analysis. 

Figure  2  shows  the  mounting  locations  of  the  UVISI  sensors  on  the  MSX  satellite.  The  satellite  is 
three-axis  stabilized  and  has  a  ring  laser  gyro /star-camera  attitude  determination  system.  All  UVISI 
sensors  are  mounted  on  the  periphery  of  the  main  spacecraft  truss  assembly  and  share  a  common 
boresight  direction.  For  target  tracking  and  attitude  control  the  entire  spacecraft  is  slewed  in  azimuth 
elevation  and  roll  by  means  of  four  momentum  wheels.  Excess  momentum  is  dumped  magnetically  by 
means  of  3-axis  control  rods.  The  spacecraft  can  be  slewed  in  azimuth  and  elevation  at  up  to  1.6Vsec 
with  accelerations  of  up  to  0.03* /sec^.  Post  processed  attitude  can  be  determined  to  1.8  arcsec  while,  orr- 
board  attitude  can  be  determined  to  0.1*.  The  on-board  attitude  is  stabilized  during  a  data  collection  event 
to  1.8  arcsec/sec. 

The  five  spectrographic  imagers  (SPIMs)  provide  two-dimensioiuil  maps  of  the  spectra  of  various 
objects  and  scenes.  Although  primary  SPIM  observations  are  the  atmospheric  dayglow  and  nigjrtglow, 
aurora,  stars,  zodiacal  light,  and  plume  trails,  it  is 


UVISI 


Fig.  1.  Block  diagram  of  UVISI  showing  the  major  subsystems  and  their  interconnections. 


SIDE  VIEW 


TOP 

VIEW 


SPIMs  (5) 


electronics  superstructure 


Fig.  2.  Mounting  locations  of  the  UVISI  sensors  on  the  MSX  spacecraft. 

All  sensors  have  a  conunon  boresig^t . 

possible  to  conduct  spectrographic  imaging  of  low  altitude  RSOs.  When  suitably  inverted,  SPIM  radiarKe 
measurements  will  reveal  atmospheric  properties  such  as  species  concentratior\s,  temperatiu’es,  and 
altitude  profiles.^'®'^  SPIM  investigations  of  auroral  radiarKe  can  provide  estimates  of  the  fluxes  and 
energies  of  precipitating  particles  that  cause  auroral  emissions  SPIM  observatioru  will  also  measure 
the  diffuse  ultraviolet  background  of  the  sky  Finally,  SPIM  observations  have  the  potential  to  provide 
composition  iitformation  on  nearby  RSOs. 

The  four  imagers  provide  both  wide-field  and  narrow-field  images,  in  selectable  wavelength 
bands,  of  the  scenes  that  the  SPIMs  observe  and  will  supply  target  tracking  information  for  use  by  the 
image  processor,  track  processor  and  satellite  control  system.  The  imagers  provide  mformation  on  the 
spatial  content  of  the  scene  to  a  resolution  of  -100  firad  or  -  20  arcsec.  Also,  the  imagers  take  one 
sitapshot  of  the  scene  every  0.5  or  0.25  seconds,  which  provides  a  much  higher  time  resolution  than  the 
SPIMs.  Thus,  tiie  imagers  can  observe  the  time-evolution  of  the  scene  and  target. 

3.1MAGEBS 

The  two  UVISI  narrow-field  imagers  share  a  common  optical  design,  shown  in  tiie  top  diagram  of 
Elg.  3.  Each  NFOV  imager  has  a  rectangular  field  of  1.6”  x  13”,  which  translates  to  a  pixel  size  of  about 
100  ^tad  in  the  256x244  format  of  the  CCD.  One  of  tire  NFOV  imagers  has  an  RbTe  photocathode  that 
resporKis  to  photons  in  the  middle  ultraviolet  (200-300  run),  which  is  the  so-called  "solar-bliiKl"  spectral 
regime.  The  other  NFOV  imager  has  an  ExtS20  (bialkali)  photocathode  that  responds  to  photons  in  the 
near  ultraviolet  arvl  visible  (300-800  nm).  Both  seirsors  use  reflecting  optics  of  a  Cassegrainian  design. 
Each  NFOV  imager  has  a  6-position  filter  wheel,  one  position  of  which  is  a  closed  or  opaque  filter.  Each 
NFOV  imager  also  has  an  effective  collecting  area  of  130  cm^. 

The  two  UVISI  wide-field  imagers  do  not  have  a  common  design.  Each  WFOV  imager  has  a 
rectangular  field  of  13.1”  x  103”,  which  translates  to  a  pixel  size  of  about  800  ^rad.  To  save  weight  and 
volume,  the  WFOV  visible  imager  (Fig.  3,  middle  diagram)  employs  a  4-element  refractive  lefts  design. 
The  WFOV  ultraviolet  imager  (Fig.  3,  bottom  diagram)  responds  to  far-ultraviolet  photons  (110-200  nm) 
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and  must  employ  fully  reflective  optics.  The  WFOV  ultraviolet  imager  uses  three-element  reflective 
optics.  As  with  the  NFOV  imagers,  each  WFOV  imager  has  a  six-element  filter  wheel  with  one  position 
opaque. 


4.  SPECTROGRAEHIC  IMAGERS 

Figure  4  shows  a  schematic  diagram  of  a  spectrographic  imager.  All  UVISl  SPIMs  have  a 
common,  off-axis  design  that  employs  three-element  reflecting  optics.  A  scan  mirror  sweeps  die  field  of 
view  through  1.0®  and  sends  light  through  a  five-position  filter/slit  mechanism.  A  collimating  mirror 
directs  the  li^t  from  the  slit  to  a  curved  grating,  which  was  mechanically-ruled.  The  spacing  of  the 
gratings  differ  for  each  SPIM  in  order  to  provide  the  desired  dispersive  characteristics.  The  grating  itself 
serves  to  focus  the  dispersed  ligj\t  on  a  two-dimensional  focal  plane.  One  dimension  of  the  focal  plane 
represents  a  spatial  dimension  of  1°,  while  the  other  dimension  represents  a  spectral  dimension  equal  to 
the  dispersion  of  the  spectrograph.  One  complete  2D  spatio-spectrogram  results  from  each  measurement 
at  each  n»irror  position.  The  five  position  slit/filter  mechanism  provides  two  slit  sizes  (1.0®x0.10®  and 
1.0®x0.05®)  and  various  blocking-filters  that  eliminate  extraneous  spectral  orders  and  long-wavelength 
contaminants. 
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Fig.  3.  Schematic  diagrams  of  the  three  types  of  UVISl  imagers 
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Figure  5  illustrates  the  corKept  of  tf\e  UVISI  spectrographic  imagmg-  The  scan  mirror  sweeps  the 
slit  throu^  1®  in  discrete  steps  of  0.10®  or  0.05®.  The  step  size  of  the  mirror  scan  can  be  programmed  to 
be  0.05®  or  0.10®,  corresponding  to  the  slit  size.  The  long  dimension  of  the  slit  is  also  1®  and  corresp>onds 
to  the  spatial  dimension  of  the  focal  plane.  The  second  dimension  of  the  focal  plane  corresponds  to  the 
spectral  or  dispersive  dimensicm  of  the  spectrograph.  The  CCD  can  be  programmed  to  sum  counts  on- 
chip  so  as  to  give  several  combinations  of  pixel  sizes  in  both  the  spatial  and  spectral  dimensions. 
Spedflcally,  the  CCD  array  can  have  272,  136,  or  68  pixels  in  the  spectral  dimension  and  40,  20,  10,  or  5 
pixels  in  tlw  spatial  dimension.  The  particular  combination  chosen  depends  on  the  telemetry  mode  of  the 
spacecraft  and  on  the  scene  being  viewed. 
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Hg.  4.  Schematic  diagram  of  a  UVISI  spectrographic  imager  (SPIM). 

5.  FOCAL  PLANE  UNITS 

Each  UVISI  sensor  has  a  focal  plaiw  unit  with  a  common  design,  as  shown  in  Hg.  6.  (The  design 
differs  only  in  photocathode  material.)  Incident  light  stimulates  photoelectrons  in  the  photocathode 
material,  and  tl^  photoelectrons  are  cascade-amplified  by  a  chevron  micro-channel  plate  (MCP).  The 
amount  of  amplification  or  "gain"  can  be  adjusted  by  changing  the  voltage  across  the  MCP.  The 
amplified  electrcxis  strike  a  green  phosphor,  which  produces  green  light  that  travels  down  a  fiber-optic 
taper  to  a  charge<oupled  device  (CCD).  The  CCD  converts  the  green  li^t  photons  into  digital  signals. 

The  CCD  provides  an  interscene  dynamic  range  of  12  bits  or  about  4  x  1(P.  Throu^  the  variation 
of  the  MCP  volta«  and  its  cycle  time,  the  image  intensifier  provides  an  additional  interscene  dynamic 
range  of  about  10^.  At  its  highest  gain,  the  focal  plane  unit  operates  in  a  photon-counting  mode  in  which 
four  pixel  output  words  ("counts")  represent  one  photo-event  at  the  photocathode.  The  CCD  outputs  a 
maximum  of  256x244  pixels  for  fire  imagers  and  272x40  pixels  for  the  spectrographic  imagers.  The 
intensifier  gain  is  selected  automatically  based  on  the  counting  statistics  of  the  raw  image;  however,  the 
gain  state  can  be  overridden  by  manual  command. 


25 


40,20,10  or  54 
Spatial  Elements 


Fig.  5.  Schematic  of  how  spectrographic  imaging  works  in  UVISI  SPIMs. 
6._WAVELENGTH  COVERAGE 

Figure  7  graphically  illustrates  the  wavelength  coverage  of  the  rune  UVISI  ser«ors.  Each 
rectangular  block  indicates  tixe  full  bandwidtii  of  one  sensor.  The  top  grouping  shows  the  imagers  and 
their  five  effective  filters;  the  bottom  grouping  ^ows  the  full  dispersion  of  each  SPIM.  The  filters  have 
been  chosen  to  select  wavebands  of  scientific  interest.  In  addition,  the  WFOV  visible  imager  (IVW)  has  a 
near-focus  filter  that  allows  examination  of  small  particles  (self  contaminants  or  co-orbiting  debris)  at 
ranges  of  a  few  meters  in  front  of  the  spacecraft.  The  NFOV  ultraviolet  imager  (lUN)  has  a  filter  at  290 
run  that  also  has  polarization  properties  that  allow  examination  of  polarized  ligjit.  Certain  SPIM  slit 
positions  have,  in  addition  to  blod^g  filters,  filters  ttrat  attenuate  unwanted  ligjrt  within  the  bandwidth 
of  the  sensor.  For  example,  SPIMl  has  a  filter  that  attenuates  unwanted  Lyman  alpha  radiation  (121.5 
nm)  and  OI  radiation  (130.4  run) ,  while  allowing  passage  of  LBH-band  photorrs  (>140  run). 
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Fig.  6.  Schematic  diagram  of  UVISI  focal  plane  unit.  Focal  plane  units  differ  only  in  material  of 

photocathode. 


7.  IMAGER  SENSmvnY 

The  purpose  of  this  discussion  is  to  address  the  suitability  of  the  UVISI  instruments  to  acquisition 
and  tracking  of  RSOs  at  both  the  visible  and  solar  blind  UV  wavelengths.  A  complete  discussion  of  the 
UVISI  instrument  sensitivity  as  calibrated  can  be  found  in  the  series  of  calibration  reports  (one  for  eadi 


Figure  8  is  a  graphic  demonstration  of  tite  inherent  utility  of  a  UV  sensor  in  detecting  objects  near 
tile  hard  earth  or  in  eailh  limb.  These  data  were  collected  by  an  imaging  instrument  flown  on  a  LEO 
satellite.  The  instrument  and  telescope  were  of  similar  design  to  tiie  UVISI  lUN.  The  object  in  the  image 
field  is  seen  in  reflected  sunlight  at  a  wavelength  of  230  nm  and  is  easily  detected  even  against  tiie 
relatively  bright  Eartii  limb  background.  The  range  to  the  object  was  >-1500  Im  at  tiie  time  that  this  image 
was  collected. 

Although  all  nine  UVISI  instruments  have  potential  utility  to  a  space  surveillance  application,  tiie 
two  instruments  of  primary  interest  to  this  community  are  the  IVN  and  the  lUN  narrow  field  of  view 
ultraviolet  and  visible  imagers.  Figures  9  and  10  display  tiie  level  of  sensitivity  required  by  the  IVN  and 
lUN  imagers  to  detect  a  Im^  plate  as  a  function  of  range  and  of  wavelength.  The  plate  is  assumed  to  be 
lambertian  and  have  a  reflectance  of  0.8. 
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The  curves  show  the  intensities  of  reflected  sunlight  convolved  with  each  wavelength  band  for 
each  of  the  various  filters.  The  horizontal  lines  represent  the  imager  sensitivity  at  that  particular 
wavelength.  Thus,  ttte  intersection  of  a  curve  with  its  respective  sensitivity  represents  the  range  expected 
for  which  a  Im^  object  can  be  detected  at  a  given  wavelength. 
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Energy  Flux  vs  Range  (IV N) 

(Solar  Spectrum  @  UVISl  tVN  Bandpass  for  1m2  Plate  with  0.8  reflectance) 
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Fig.  9.  Intensity  of  reflected  sunlight  vs.  Range  of  a  1  plate  with  an  albedo  of  0.8.  \  is  binned  to  the 
IW  filter  bands  and  the  respective  sensitivity  is  shown  by  horizontal  lines . 
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Fig.  10.  Intensity  of  reflected  sunlight  vs.  Range  of  a  1  plate  with  an  albedo  of  0.8.  X  is  binned  to  the 
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Energy  Flux  vs  Range  (lUN) 


( Solar  Spectrum  Q  UVISl  lUN  Bandpass  for  1m2  Plata  with  0.8  reflectance) 
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Figure  9  shows  that  the  IVN  Open  (300-900  nm),  2  (350-440  nm)  and  3  (470-640  nm)  filter 
positions  can  detect  a  Im^  plate  which  lies  beyond  geosynchronous  orbit.  In  the  filter  1  (305-315  run  ) 
position  IVN  is  limited  to  a  range  of  -  6000  km  for  detection  of  the  same  Im^  plate. 

Figure  10  shows  that  the  detection  range  for  the  filter  positions  of  the  lUN  varies  widi 
wavelength  from  -  20(XX)  km  for  filter  1  (180-300  run),  about  150(X)  km  for  filter  3  (260-300)  and  falls  off  to 
-  4000  km  for  filter  2  (200-230  nm). 

8.imageprcx:essor 

The  image  processor  of  UVISI  makes  the  instrument  unique  among  those  flown  to  date.  The 
image  processing  system  can  acquire  and  isolate  likely  targets  in  an  image  FOV  and  communicate  their 
positions  to  the  MSX  flight  processor^  which  in  turn  points  the  satellite  in  a  "closed-loop"  fashion. 
"Target"  can  refer  to  either  a  point  source  such  as  a  star  or  other  satellite  or  to  a^t  extended  source  such  as 
an  auroral  or  cloud-top  feature.  The  image  processor  operates  in  real  time  and  can  automatic2dly  select 
the  operating  mode  of  the  UVISI  imager  based  on  its  results. 

The  image  processor  accepts  raw  image  data  (a  256x244  pixel  array  of  counts)  from  any  orte  of  the 
UVISI  imagers.  Processor  software  performs  initialization  and  tracking  fimctions  fiiat  include  filtering, 
smoothing,  thresholding,  and  centroiding.  The  image  processor  seeks  targets  based  on  an  a  priori  target 
description  file  containing  weights  for  numerous  target  features  such  as  size,  shape,  brightness,  and 
location.  The  image  processor  then  traitsforms  target  locations  from  UVISI  pixel  coordinates  to  spacecraft 
coordmates  and  passes  titem  to  the  MSX  flight  processor,  which  performs  Kalman  filtering  of  other 
targeting  inputs  to  select  "true"  targets  for  observation. 

The  image  processor  offers  considerable  flexibility  in  its  operation.  The  nomitud  processor 
program  determines  likely  targets  based  on  a  feature  file  contained  in  its  memory.  The  software 
maximizes  target  likelihood  by  weighting  the  features  in  tibe  file.  One  can  change  target  selection  by 
changing  the  weights  in  the  file.  Thus,  a  star  target  would  be  weighted  to  optimize  bri^tness  and 
pointedness,  while  tite  active  segment  of  an  aurora  arc  would  be  weighted  to  optimize  elongation  of 
shape  and  orientation.  Because  more  tiian  a  single  target  can  appear  in  one  frame  of  data,  tite  image 
processor  can  accommodate  up  to  23  different  targets  at  cmce. 

UVISI's  image  processor  will  initially  carry  a  default  set  of  target  description  vectors, 
characterizing  a  representative  set  of  targets  by  such  properties  as  average  bri^tness  and  eccentricity. 
For  some  of  the  targets  and  phenomena  to  be  observed,  these  descriptions  are  not  much  more  thw 
educated  guesses.  During  the  early  operations  phase,  analysis  of  UVISI's  image  data  will  provide 
guidance  in  developing  better  description  vectors.  These  will  be  tested  with  tire  tracking  processor  to 
improve  overall  system  performaiKe. 

9.  TRACKING  fe  TRACK  PROCESSOR 

The  fundamental  design  driver  of  the  MSX  tracking  system  are  the  various  observatioi\s  to  be 
made  by  the  spacecraft  sensors.  Of  tiiese  the  most  hig^y  valued  observations  involve  the  launch  of 
missiles  and  sounding  rockets  dedicated  to  the  MSX  mission.  Other  high  priority  targets  irKlude 
cooperative  missile  launches.  Resident  Space  Objects,  and  reference  objects  to  be  deployed  from  canisters 
onboard  MSX.  Background  observations  of  die  Eardt  limb  and  calibration  observations  of  stars  will  also 
be  conducted.  One  additional  observation  which  has  elements  of  both  background  observations  and 
target  tracking  is  tiie  autorwmous  tracking  of  auroral  bright-spots  as  tiiey  develop  and  move  through  the 
polar  limb. 
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The  fundamoital  requirement  that  flows  down  to  the  tracking  system  from  this  observation  set  is 
therefore,  to  provide  a  continuous  "best"  estimate  of  target  state  at  a  rate  compatible  with  the  attitude 
control  system's  requirements,  whether  or  not  sensor  data  is  provided  to  or  used  by  the  track  processor. 

In  order  to  accomplish  this,  a  loosely  coupled,  distributed  processing  approach  to  tracking 
architecture  was  chosen,  partly  because  no  single  space-qualified  processor  was  foimd  to  be  capable  of 
the  total  computational  load,  aind  partly  to  reduce  developmental  risk.  Figure  11  is  a  block  diagram  diat 
schematically  represents  the  pointing  control  system.  The  shaded  boxes  focus  on  tracking  and 
experiment  control. 


Figure  11  Spacecraft  Pointing  Control  Architecture 


Hgure  12  Tracking  system  Hardware  Context 


Figure  12  shows  the  hardware  associated  with  tracking  and  the  interfaces  between  the 
components.  In  this  figure,  each  of  the  shaded  boxes  represents  a  spacecraft  subsystem,  and  each 
arrowhead  corresponds  to  a  hardware  interface.  The  boxes  shown  for  the  sensors  are  dtose  parts  of  the 
sertsor  systems  widt  direct  interfaces  to  the  tracking  processor;  UVISrs  Image  Processor,  SPIRITS  On¬ 
board  Signal  and  Data  Processor,  and  SBV's  Experiment  Controller. 
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Tracking  and  Attitude  Processors 

The  tracking  processor  and  the  attitude  processor  have  identical  central  processing  boards,  built 
by  APL  around  a  Performance  Semiconductor  1750A  PACE  chipset.  The  attitude  processor  has  a  dual- 
CPU-board  system  witii  a  serial  port  between  them.  Each  of  the  CPU  boards  has  256K  words  of  error- 
detecting-and*correcting  RAM,  a  2  K  word  boot  PROM,  a  memory  management  unit,  arxi  a  processor 
interface  circuit.  An  expansion  board  provides  256  K  words  of  EEPROM  used  for  storage  of  code, 
constants  aiKl  data  to  be  retained  between  power  cycles.  Also  on  the  expansion  board  are  an  additiorud 
256  K  words  of  RAM  and  an  interrupt  expansion  chip.  The  CPU  is  rated  at  1.3  MIPS  (DAIS  mix  with 
floating  point)  at  20  MHz  clock  rate.  In  addition  to  titese  boards,  a  variety  of  custom  interface  boards  are 
in  each  computer  chassis.  Most  of  the  custom  interfaces  are  based  on  FIFO  buffer  circuits,  with  interrupts 
provided  to  the  CPUs  at  the  completion  of  a  predetermined  number  of  words  for  each  interface.  The 
tracking  and  attitude  processors  are  fully  redundant,  as  are  all  spacecraft  life-critical  subsystems. 

Data  Flow 

Refer  to  figure  12  to  follow  the  data  flow  through  the  tracking  system.  The  attitude  processor 
passes  a  time-tagged  message  to  the  tracking  processor  twice  per  second.  The  message  contains  the 
estimated  attitude  quaternion,  attitude  body  rates  and  accelerations,  the  spacecraft  orbital  position  and 
rate,  and  status  information  about  the  attitude  processor.  The  tracking  processor  bases  its  response  on 
control  parameters  received  from  the  commarvl  system  that  describe  the  nominal  trajectory.  It  returns  to 
the  attitude  processor  a  time-tagged  message  witii  a  desired  attitude  quaternion,  rate  and  acceleration, 
and  status  information  about  tire  tracking  processor.  The  rate  for  this  message  is  also  2  Hz. 

The  desired  attitude  state  may  be  based  entirely  on  tire  nortrirral  trajectory,  or  it  may  be  based  on 
setrsor  observatiorrs.  To  support  the  sensor  observatiorrs,  tire  tracking  processor  provides  at  2  Hz  a  single 
rtressage,  broadcast  by  hardware  to  UVISI,  SBV  and  OSDP.  The  message  contairrs  a  target  state  estimate 
in  Earth-centered  inertial  coordinates  and  again  in  UVISI-referetKed  coordirrates  (ffactional  pixels  of  the 
currently  selected  imager).  The  spacecraft  attitude  and  orbit  data  from  the  attitude  processor  are  also  sent 
in  the  message,  along  with  some  information  about  the  expected  observatioiral  conditions  of  the  target, 
and  status  iirformation  wr  the  tracking  processor.  The  nressages  from  the  sensors  to  the  tracking  processor 
are  discussed  below  in  the  sectirm  "Closed-Loop  Tracking.” 

Specification  and  Control  of  Data  Collection 

A  data  collecticm  event  may  comprise  any  number  of  tracking  events,  witii  three  perhaps  typical 
(calibrate,  observe  target,  calibrate).  Each  of  ^ese  tracking  events  is  described  separately,  and  the 
parameters  passed  to  the  tracking  processor  itKlude  the  start  and  stop  times  for  each.  While  oi>e  event  is 
executing,  messages  stored  in  ^e  command  processor  memory  may  be  transferred  to  tire  tracking 
processor,  with  start  time  for  the  new  event  preceding  stop  time  for  the  current  one.  In  ttris  way  an 
indetinite  chaining  of  tracking  events  may  be  accomplished,  with  smooth  transition  ffom  one  to  the  next, 
and  no  return  to  parked  mode  in  between.  Blending  of  trajectories  from  the  end  of  one  event  to  the 
begiiming  of  the  next  is  the  responsibility  of  operations  team  planners,  not  an  ort-board  hmcticm. 

Tracking  event  control  parameters  are  broken  into  three  categories:  those  that  are  usually 
ccmstant,  like  sensor  alignment  matrices;  those  that  define  the  next  event,  such  as  the  iwminal  trajectory; 
and  those  that  are  intended  to  change  the  way  tiie  current  event  is  being  executed  such  as  a  revision  to 
the  choice  of  sensors  to  be  used  for  tracking  the  currott  target.  Any  of  ^se  parameters  may  be  sent  in 
real  time,  but  generally  the  parameters  will  be  stored  in  the  command  system  for  delivery  to  foe  tracking 
processor  at  predetermined  times. 

When  control  parameters  are  sent,  they  may  be  stored  in  EEPROM  to  become  the  new  defaults. 
Then  whenever  foe  tracking  processor  is  rebooted,  it  will  load  those  control  parameters.  Alternatively, 
the  parameters  may  be  stored  only  in  RAM,  and  their  values  will  be  lost  when  foe  processor  is  turned  off 
after  the  event.  If  an  experiment's  control  parameters  are  stored  in  EEPROM,  foe  processor  may  be 
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turned  on  by  stored  command  even  if  it  out  of  view  of  any  ground  station,  and  the  experiment  can  be  run 
without  any  real-time  contact  with  the  tracking  processor. 

Tracking  Modes 

Many  data  collection  events  do  not  require,  or  could  not  use,  sensor  data  to  improve  pointing 
control.  All  that  is  needed  is  a  specification  of  the  desired  maneuvering  as  a  function  of  time,  i.e.  open 
loop.  The  sensors  then  send  their  data  to  the  tape  recorder,  but  iv>  use  is  made  of  it  on  board.  While  this 
kind  of  tracking  event  is  simpler  than  a  closed -loop  track,  there  are  still  many  optiorts.  Furthermore,  any 
closed-loop  tracking  event  requires  a  nominal  trajectory  as  the  basis  for  the  track  and  to  specify  where  the 
target  is  to  be  acquired  14  . 

When  dwre  is  a  sufficiently  detectable,  compact  target,  especially  orw  whose  dynamics  are  kitown 
fairly  well,  closed-loop  tracking  can  be  used  to  advantage.  LFVISI  and  the  tracking  processor  cooperate  in 
closed  loop  target  tracking.  The  UVISl  image  processor  examines  a  single  hame  from  one  of  its  four 
imagers,  and  sends  to  the  tracking  processor  a  rank-ordered  list  of  the  objects  it  detects.  The  frame  rate, 
and  the  message  rate,  is  2  Hz.  The  objects  are  ranked  according  to  their  goodr>ess-of-fit  to  a  feature  vector 
uplinked  to  dte  image  processor. 

The  list  of  potential  targets  from  the  image  processor  may  be  used  in  either  of  two  ways,  as 
specified  by  experiment  control  parameters.  The  tracking  processor  may  be  instructed  simply  to  use  the 
object  rartked  l^ghest  by  the  image  processor  that  is  within  a  reasonable  distance  of  the  estimated  target 
location.  On  the  other  hand,  the  tracking  processor  can  use  its  velocity  filtering  algoridun  to  choose  dtat 
object  in  die  list  that  most  closely  matches  the  velocity  profile  of  the  target,  using  the  image  processor's 
rankings  to  break  ties.  This  is  us^  for  star  removal. 

The  t}rpe  of  closed-loop  tracking  to  employ  depends  primarily  on  the  a  priori  knowledge  of  target 
dynamics.  If  little  or  nothing  is  known  about  the  forces  governing  target  motion  in  advance,  then  it  is 
appropriate  to  use  die  alpha-beta  tracker.  This  algorithm  predicts  the  future  IcKation  and  velcKity  of  the 
target  by  assigning  weight  alpha  to  the  position,  beta  to  the  velocity,  and  linearly  propagating  the  state  in 
these  two  variables.  This  algorithm  is  sometimes  called  a  ’data  follower,’  because  it  accumulates  no 
confidence  in  its  predictions  based  on  diminishing  error,  and  uses  no  underlying  model  of  die  target's 
dynamics.  It  has  minimal  risk  of  diverging  from  d«  target  unless  the  sensors  stop  providing  reasonable 
data,  or  target  dynamics  exceed  spacecraft  slewing  capabilities,  but  it  has  very  liMe  ability  to  follow  die 
target  througjh  extended  data  dropouts. 

A  second  method,  the  extended  Kalman  filter,  bases  its  state  propagation  on  underlying 
equations  of  motion  for  the  target.  As  sensor  data  is  received,  the  model's  state  vector  is  correct^  to 
minimize  differences  between  it  and  the  observed  state.  When  the  predicted  state  agrees  closely  with  the 
sensor  data,  confidence  is  increased  that  the  state  vector  convergence  is  good.  With  ii^udidous  setting  of 
the  filter  gains,  convergence  time  may  be  artificially  increased,  or,  in  the  catastrophic  case,  sensor  data 
may  be  r^ected  in  fovor  of  erroneous  model  predictions,  and  ttie  filter  may  actually  diverge. 

With  the  proper  settings,  however,  an  extended  Kalman  filter  can  smootii  sensor  noise  effects, 
minimize  tiie  effects  of  "wild  points,"  and  propagate  target  state  accurately  throu^  times  of  sensor  data 
dropout.  Particularly  when  the  target  dynamics  are  very  well  known,  as  in  ballistic  motion,  the  Kalman 
filter  can  effectively  correct  wrong  initial  assumptions  about  the  target's  positim  and  velocity,  converging 
to  a  state  vector  that  minimizes  sensor  data  residuals. 

If  closed  loop  tracking  is  to  be  performed,  one  or  more  sensors  must  be  selected.  The  tracking 
processor  may  be  commanded  to  use  one  sensor  until  and  unless  it  loses  tiie  target,  and  tiien  switch  to  the 
next  priority  sensor,  or  it  may  be  told  to  give  two  or  all  three  of  the  sensors  the  same  priority.  If  two  or 
more  sensors  are  selected  for  simultaneous  use,  s«isor  fusion  is  specified.  When  incorporating  sensor 
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data  using  either  the  Kalman  filter  or  tfw  alpha-beta  tracker,  it  is  effectively  weighted  by  its  frequency  of 
arrival,  in  the  sense  that  each  data  point  from  a  source  is  given  equal  weight,  without  regard  to  its  rate. 

Loss  of  track  by  all  selected  sensors  leaves  the  tracking  processor  with  two  choices:  continue  to 
propagate  the  last  predicted  target  state,  or  revert  to  the  nominal  trajectory.  The  decision  is  made  in 
advance  by  spacecraft  operators.  Note  that  the  tracking  processor  may  be  operated  in  either  alpha-beta 
track  or  in  Kalman  filter  mode  while  the  spacecraft  is  pointed  inertiaJly.  In  this  mode  of  operation  target 
track  files  can  be  collected  on  objects  which  pass  through  a  sensor's  field  of  view  and  at  each  loss  of  track 
a  new  target  will  be  selected. 

Sensor  Fusion 

Data  from  very  different  sensors  may  be  combined  into  a  single  target  state  estimate.  There  is 
much  to  learn  about  the  process.  When  incorporating  the  sensor  data  into  the  prediction  update,  a  model 
of  sensor  noise  is  applied.  By  increasing  the  expected-noise  values  for  a  particular  sensor,  the  influence  of 
the  data  may  be  decreased.  Thus,  by  adjusting  the  noise  models,  same-priority  sensors  with  different 
message  frequencies  can  be  equalized.  For  targets  diat  have  distinctive  signatures  in  one  spectral  band, 
intuition  says  that  the  weight  of  corresponding  sensor  data  should  be  increased.  The  MSX  tracking 
system  can  ^ow  us  to  examine  this  aivi  other  issues  in  multispectral  sensor  fusion. 

Another  way  to  use  multiple  sensors  for  one  tracking  event  is  to  switch  from  one  to  another  at  a 
predetermined  time.  Some  phenomena  may  have  definable  points  in  their  trajectories  at  which  such  a 
change  of  sensors  is  called  for.  For  example,  the  accuracy  provided  by  the  beacon  receiver  au%  not  as 
great  as  for  UVISI.  It  may  make  sense  therefore  to  start  tracking  a  cooperating  test  v^cle  with  fite 
beacon  receiver,  aird  then  switch  over  to  UVISI  at  some  point  in  the  vdticlc’s  flight.  For  such  highly 
valued  targets,  however,  it  is  certaiidy  advisable  to  quantify  the  benefits  of  this  approach  first 

Simulation 

For  testing  of  the  tracking  system  software  at  higgler  levels  of  integration,  a  tesfibed  simulator^ 
was  developed.  The  simulator  provides  full  electrical-level  compatibility  at  all  tracking  system  interfaces, 
artd  varying  degrees  of  fidelity  in  the  data  exchanged  across  them.  It  allows  different  degrees  of 
hardware-in-the-loop  simulation,  including  the  tracking  processor,  attitude  processor,  beacon  receiver, 
and  UVISI  image  processor.  The  simulator  has  supported  test  generation,  data  capture  aiwl  analysis 
throughout  the  development  effort,  and  will  be  used  after  launch  for  system  characterization  and 
software  mainterumce. 

10.  SUMMARY 

The  preceding  paper  was  meant  to  provide  file  reader  wifii  some  insight  into  the  capabilities  of 
the  UVISI  Instrument  coupled  with  its  Image  and  Track  processors.  It  is  clear  fiiat  target  selection,  track 
file  computations  and  trajectory  determination  can  be  accomplished  in  a  variety  of  ways  and  at  various 
different  wavelengths.  The  UVISI  Narrow  Field  Imagers  are  best  suited  to  collect  track  data  on  RSOs  at 
ranges  up  to  and  beyond  geosyndironous.  The  overall  dynamic  range  of  ~10^  and  the  solar  blind  UV 
wavelength  range  allows  UVISI  to  view  RSOs  and  other  dim  objects  in  deep  space  as  well  as  objects  in 
the  earth  limb  or  against  a  sunlit  earth  background.  When  coupled  with  fiie  MSX  tracking  and  control 
system,  UVISI  has  the  potential  to  acquire  and  track  RSOs  providing  necessary  track  file  information  and 
long  duration  exposures.  Spectral  information  from  file  SPIMs  could  thus  be  collected  and  contribute  to 
multispectral  SOI  by  file  identification  of  specific  reflectaxKe  spectra. 


34 


11.  ACKNOWLEDGMENTS 

The  authors  would  like  to  thank  the  large  number  of  supporting  engineers  and  scientists  who 
contributed  to  the  developntent  and  construction  of  UVISI  aitd  the  Image  and  Track  Processors.  The 
UVISI  instrument  and  the  MSX  mission  receives  funding  through  the  Ballistic  Missile  Defense  Office, 
formerly  the  Strategic  Defense  Initiative  Organization. 


9.  REFERENCES 

1.  C.-I.  Meng  and  R  .E.  Huffman,  "Ultraviolet  imaging  from  space  of  the  =iuroral  under  full 
sunlight,"  Geophys.  Res.  Lett.,  11, 315-318, 1984. 

2.  F.  W.  Schenkel  and  B.S.  Orgorzalek,  "Auroral  images  from  space:  Imagery,  sj)ectroscopy,  and 
photometry,"  APL  Tech.  Dig.  8, 309-317, 1987. 

3.  A.F.  Davidson,  K.S.  Long,  S.T.  Ehirrance,  W.P.  Blair,  C.W.  Bowers,  S.J.  Conard,  P.D.  Feldman, 
H.C.  Ferguson,  G.H.  Fountain,  R.A.  Kimble,  G.A.  Kriss,  H.W.  Moos,  and  K.A.  Potocki,  "The  Hopkins 
Ultraviolet  Telescope:  Performance  and  calibration  during  the  Astro-1  Mission,"  Astnphys. ).  392, 264-271, 
1992. 

4.  M.  D.  Grifhn  and  M.  j.  Rendirte,  "Delta  180/Vector  Sum:  the  first  powered  space  intercept,"  AIAA 
26th  Aerospace  Sciences  Meeting,  Reno,  NV,  January  1988. 

5.  B.  B.  Holland,  J.  C.  Ray,  and  A.  J.  Green,  "Delta  181,  the  Thrusted  Vector  Mission,"  AIAA  27th 
Aerospace  Sciences  Meeting,  Reno,  NV,  }anua.y  1989. 

6.  P.  F.  Bythrow,  J.  F.  Carbary  and  R.  E.  Gold,  Delta  Star  Target  Observatiorts  at  Wavdengths  <  1pm, 
Proceedings  of  the  Short  Wavelength  Phenomenology  Conference  JHU/APL  June  1990. 

7.  S.  C.  Solomon,  P.  B.  Hays,  and  V.  J.  Abreu,  "Tomographic  inversion  of  satellite  photometry," 
Appl.  Opt.,  19, 3409-3418, 1985. 

8.  R.  R.  Meier,  "Ultraviolet  spectroscopy  and  remote  sensing  of  the  upper  atmosphere,"  Space  Sci. 
Reo.,  58, 1-185, 1991. 

9.  D.  J,  Strickland,  R.  Link,  and  L.  J.  Paxton,  "Far-UV  remote  sertsing  of  thermospheric  composition 
and  the  solar  EUV  flux,"  in  Ultraviolet  Teclmology  FV,  R.  E.  Huffman,  ed.,  Proc.  Soc.  Photo-Opt.  Instr.  Eng., 
1764, 117-131, 1992. 

10.  D.  J.  Strickland,  J.  R.  Jasperse,  and  J.  A.  Whalen,  "DependeiKe  of  auroral  FUV  enussions  on  the 
irKident  electron  spectrum  and  neutral  atmosphere,"  /.  Geophys.  Res.,  88, 8051-8062, 1983. 

11.  R.  C.  Henry,  J.  Murthy,  M.  Allen,  M.  Corbin,  and  L.  J.  Paxton,  "Spectroscopy  and  imaging  of  the 
cosmic  diffuse  UV  background  radiation,"  in  Ultraviolet  Technology  IV,  R.  E.  Huffman,  ed.,  Proc.  Soc. 
Photo-Opt.  Instr.  Eng.,  1764, 61-71, 1992. 

12.  J.  F.  Carbary  The  UVISI  Calibration  Report,  JHU/APL  Publication  1994. 

13.  L.  Lee  Pryor  and  Daniel  S.  Wilstm,  MSX  Flight  Software  Test  Bed 

14.  R.  L.  Waddell  Jr.,  P.  K .  Murphy  and  G.  A.  Heyler,  Image  and  Track  Processing  in  Space,  Part  1, 
Proc.  of  the  AIAA  conference  on  Computing  in  Aerospace,  pp.  576-585  Oct.  1993. 


% 


35 


The  SBV  Sensor  Emulator;  A  High  Fidelity  Focal  Plane  Simulator 
Thomas  P.  Qpar  ( MTT  Lincoln  Laboratory ) 

1.0  Introduction 

The  Space  Based  Visible  (SBV)  sensor  on  the  Ballistic  Missile  Defense  Organization 
Midcourse  Space  Experiment  (MSX)  satellite  collects  and  processes  image  data  from  a  staring, 
CCD  array.  As  an  aid  to  visualization  of  the  images,  provide  test  cases  for  software  and  algorithm 
development,  and  to  serve  as  the  basis  for  phenomenology  modelling,  the  SBV  Sensor  Emulator 
has  been  developed.  The  purpose  of  the  SBV  Sensor  Emulator  is  to  provide  the  SBV  user 
community  with  high  fidelity,  simulated  images.  The  SBV  program  is  already  rich  in  field 
measurement  fix)m  an  SBV-like  sensor  which  provides  one  set  of  opportunities  to  test  software. 
The  existing  field  test  data  is  useful  for  assessing  the  functionality  of  existing  software,  but  does 
not  provide  a  complete  software  assessment  due  to  the  differences  between  Ae  existing  collection 
systems  and  SBV.  There  are  six  primary  differences  which  are  specifically  addressed  by  the  SBV 
Sensor  Emulator,  which  are: 

1 . )  The  spatial  structure  of  the  multiplicative  fixed  pattern  due  to  optical  and  CCD 
effects  is  not  present  in  the  field  data. 

2. )  The  SBV  sensor  has  a  large  degree  of  optical  distortion  which  is  absent  in  the 

current  field  data. 

3 .  )  The  noise  performance  is  a  sensor  characteristic  which  is  unique  to  each  individual 

sensor. 

4.  )  The  spectral  response  of  the  sensor  differs  from  previous  measurements,  primarily 

as  a  result  of  the  SBV  mirror  spectral  reflectance. 

5 . )  The  Optical  Point  Spread  Function  (PSF)  is  sharper  than  most  of  the  field  test  data 
and  varies  across  the  focal  plane. 

6. )  The  SBV  Sensor  is  space-based  and  opens  a  variety  of  background  phenomenology 

and  encounter  RSO  encounter  geometries  which  are  not  accessible  from  the  earth.  For  example, 
the  dominant  background  contribution  for  low  tangent  altitude  observations  is  fix>m  Non-Reject^ 
Earth  Radiance  (NRER)  which  is  absent  from  existing  data. 

The  SBV  Sensor  Emulator  provides  the  means  to  alleviate  the  deficiency  in  the  available 
field  test  data  by  utilizing  the  extensive  set  of  ground  calibration  and  characterization  data  to 
provide  the  means  to  create  high  fidelity,  SBV-like  images  for  the  MSX/SBV  softv/are  users  and 
developers. 

2.0  SBV  Sensor  Emulator  Overview 

Since  this  emulator  is  being  desired  specifically  for  the  SBV  sensw  an  important  caveat 
should  be  noted.  The  SBV  sensor  is  designed  primarily  for  above  -  the  -  horizon  (ATH) 
observations  only.  This  design  implies  that  the  emulator  need  only  include  ATH  Irackground 
contributions. 

The  currait  implementation  of  the  SBV  Sensor  Emulator  includes  four  major  components: 

1 .  Cover  Closed  (Dark)  Frame  Simulation 

2 .  Extended  Source  Background  Modelling 

3 .  Stationary  and  non-Stationary  Point  Target  Simulatior 

4.  Analog  -  to  -  Digital  Conversion 

The  output  from  the  emulator  is  consistent  with  existing  SBV  data  processing  software.  It 
produces  files  which  are  in  the  format  used  by  the  SBV  CONVERT  software,  the  SPOCC 
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Processing  Pipeline,  and  the  SBV  Performance  Assessment  Team  software. 

This  paper  reviews  the  four  major  components  of  the  emulator.  The  purpose  of  this  paper 
is  to  provide  an  overview  of  the  features  and  capability  of  the  emulator.  The  maAematical 
formulation  of  the  various  elements  is  presented  in  the  following  sections.  Section  3  provides 
details  on  the  Cover  Closed  or  Dark  Frame  Simulation.  Section  4  is  a  review  of  the  Extended 
Source  Modelling.  The  techniques  to  generate  point  sources  on  the  focal  plane  are  addressed  in 
Section  5.  The  simple  Analog  -  to  -  digital  conversion  routine  is  described  in  section  6.  Section  7 
presents  the  summary  and  conclusions. 

3.0  Cover  Closed  (Dark)  Frame  Simulation 

The  emulator  begins  by  creating  the  contributions  from  the  Zero  Correction,  the  Dark 
Current,  and  the  noise.  The  Zero  Correction  is  the  parameter  which  is  used  to  set  the  DC  bias  level 
in  the  camera.  The  DC  bias  level  at  the  output  is  termed  the  Zsum  and  corresponds  to  the  mean 
signal  level  in  a  pixel  with  no  photons  applied.  The  Zsum  is  a  fiinction  of  the  Zero  Correction  and 
it’s  value  depends  on  the  specific  camera  chosen.  (The  SBV  sensor  has  two  independent, 
redundant  camera  Each  is  supported  by  the  emulator.)  TTie  program  begins  by  setting  the  signal 
level  in  each  fiame  to  the  Zsum.  This  is  done  as  a  noiseless  process  and  the  actual  operation  is 
depicted  by  Equation  1. 

s"  =  Zsum  (ZC, Camera) 

ij  1  ) 

where 


SijH  =  Signal  in  pixel  (ij)  for  frame  numbCT  n 

ZSum  =  Signal  level  for  zero  photons  as  a  function  of  the  Zero  Correction  and  the  C^era 
The  relationship  between  Zsum  and  Zero  Correction  is  linear  and  can  be  written  as 


Zsum  =  Zsum  (ZC  -  ZC^j.) 


2.) 


Table  1  lists  the  set  of  these  parameters  which  was  derived  fiom  the  SBV  characterization.  One  set 
of  parameters  is  supplied  for  each  of  the  SBV  camera  halves. 

Table  1 

Model  Parameters  for  the  Zsum  Dependence  on  Zero  Correction 


Camara 

2ium0 

ZQrf 

A 

7.93 

120.7 

B 

7.84 

121.6 

After  the  2^um  value  has  been  added  to  the  pixel,  then  the  dark  current  signal  is  added.  The  dark 
current  signal  depends  on  the  exposure  time,  while  the  dark  current  is  a  function  of  the  focal  plane 
temperature.  The  dark  current  signal  also  exhibits  a  fixed  pattern  behavior  in  which  case  the 
average  dark  current  varies  pixel  -  by  -  pixel.  The  focal  plane  distribution  of  dark  current  is 
rqrresented  by  a  two-dimensional  420  by  420  array  for  each  of  the  four  CCD  arrays.  The 
d^ndence  on  focal  plane  temperature  is  set  by  an  empirical  model.  The  dark  current  model  is 
explicitly  defined  in  Equation  3. 
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3.) 


where 


D^. 

y 


C^/T  -  T*^ 
Sj^^FPA  ^FPA 

DM.,  e 
y 


Dyn 

DMijk 

Cij" 

Tfpa“ 

Tfpa 

TEj5>^*f 

Tfixp 


=  Dark  Current  signal  in  pixel  (ij)  in  frame  n 
=  Dark  Current  in  pixel  (iJ)  on  CCD  #  k 
=  Temperature  correction  parameter  for  CCD#k 
and  pixel  (ij) 

=  Reference  focal  plane  temperature  ( ®K ) 

=  Focal  Plane  Temperature  ( °K )) 

=  Reference  focal  plane  exposure  time  (s) 

=  Focal  plane  exposure  time  (s) 


The  temperature  correction  parameter  is  set  to  a  constant  across  all  CCDs  and  pixels  in  this 
version. 

The  final  contributor  to  the  cover  closed  signal  is  the  noise.  The  dark  noise  has  two  contributions 
one  is  from  the  Poisson  Noise  due  to  the  dvk  current  signal  and  the  other  is  the  sum  of  all  other 
electronic  noises  in  the  system.  In  this  version  of  the  emulator  the  dark  noise  is  modeled  by  a  zero- 
mean  Gaussian  distribution.  The  variance  of  the  nmse  is  assumed  to  be  a  constant  value  and  is  seat 
to  ims  value  equaling  5.0  electrons  at  high  gain.  The  model  is  implemented  according  to  equation 
4. 


sD  =  s;  + 


4.) 


where  <^j  is  the  standard  deviation  of  the  noise  in  the  pixel  (i  j)  and  r(p,a)  represents  a  Gaussian 

random  variable  with  mean  p.  and  standard  deviation  a.  The  noise  value  is  detramined  by  the  root 
sum  of  squared  values  of  the  Dark  current  Poisson  Noise  and  the  system  noise: 


where 


2  2 

o..  =  <5^ 
y  sys 


+  D../g 

ij  ® 


5.) 


(Tsys2  =  system  noise  level  in  digital  numbers 
Dij  =  dark  current  signal  in  pixel  (i  j) 

g  =  system  transfer  function  or  gain  (photo-electrons  /  digital  number ) 

The  relationship  between  dark  ciurent  signal  and  the  noise  is  derived  from  the  assumption  that  the 
process  is  defined  by  Poisson  statistics. 

4.0  Extended  Source  Background  Modelling 

The  SBV  Sensor  Emulator  currently  includes  three  primary  background  constituents; 
diffuse  galactic  backgrounds,  Out-of-Field~of-View  Radiance  (OFVR)  backgrounds,  and  discrete 
celestial  point  sources.  This  section  details  the  modelling  of  the  first  two  of  Aese  contributions. 
The  last  is  described  in  Section  5. 
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4.1  Diffuse  Galactic  Background 

The  signal  in  a  pixel  due  to  a  extended  source  background  feature  that  fills  the  pixel  is  written  in  the 
following  form 


B.. 

y 


.  Eff  Optics  , ,  - 
A..  qe..T  ^  N«  10 

ij  ^  ij  0 


Tp 

ij  Exp 


6.) 


where 

Bij  =  signal  in  the  pixel  (iJ)  due  to  the  extended  background  source 

AEffy  =  effective  aperture  area  for  the  pixel  (ij) 

qeij  =  solar  weighted  quantum  efficiency  for  the  pixel  (i  j) 

xQptics  =  optical  throughput  of  the  sensor 

No  =  zero  magnitude  flux 

Bvm  =  background  visual  magnitude  per  aicsec2 

(Xij2  =  pixel  steradiancy  ( arcsec2) 

Tpyp  =  frame  exposure  time 

The  celestial  background  in  the  SB  V  Emulator  assumes  a  constant  background  visual 
magnimde  across  the  foc^  plane.  This  diffuse  background  component  is  primarily  a  function  of 
the  galactic  latimde.  It  is  due  to  the  signal  fiom  unresolved  stars  and  extended  features.  Table! 
provides  an  estimate  of  the  diffuse  gakctic  background  for  the  SBV  sensor. 

Table! 

Diffuse  Galactic  Background  Magnitude  vs.  Galactic  Latitude 


Galactic 

Bvm 

Latitude 

(mag  /arcsec2) 

(degree) 

±5 

!!.! 

±30 

!4.0 

±60 

!5.0 

The  background  level  for  which  the  Poisson  noise  in  the  background  equals  the  sensor  noise  is 
!!.!  magnitudes  /arcsec^. 

The  product  AEffy  ctjj^  is  a  function  of  the  design  of  the  sensor.  The  SBV  sensor  has  structure 
introduced  fimn  each  of  these  factors.  The  structure  in  the  multiplicative  map  is  primarily  due  to 
two  aspects  of  the  off-axis  design  of  the  SBV  sensor.  First,  the  effective  aperture  area  of  the 
telescope  varies  across  the  foc^  plane.  This  “pupil  distortion”  produces  a  radial  structure  in  the 
map  centered  on  the  on-axis  point  Secondly,  the  f/number  of  the  telescope  varies  fi’om  point  -  to  - 

point  The  variation  in  f/  number  causes  the  pixel  steradiancy  (ctij^)  to  be  smaller  in  one  portion  of 
the  focal  plane  with  respect  to  another.  Thus  when  illuminated  by  a  uniform  source  there  is  a 
gradient  in  the  signal  on  the  focal  plane. 
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4.2  OFVR  Backgrounds 


The  primary  component  of  the  background  for  low  earth  satellite  observations  may  be  due 
to  OFVR  contributions,  most  notably  the  N^R.  One  ccnnponent  of  the  NRER  is  directly  related 
to  mirror  contamination.  Quantification  of  the  contamination,  via  measiuements  of  the  Bi¬ 
directional  Reflectance  Distribution  Function  (BRDF)  of  the  SBV  mirrors,  has  been  made  over  a 
long  period.  The  BRDF  data  can  be  used  to  compute  the  expected  NRER  background  as  function 
of  tangent  altitude  for  a  fully  illuminated  earth.  In  the  SBV  Sensor  Emulator  the  NRER  is  nrodeled 
by  blowing  the  quantity  Bvm  in  Equation  6.)  to  be  a  function  of  the  tangent  altitude.  Table  3  is 
listing  of  the  expected  NRER  background  visual  magnitude  at  various  tangent  altitudes,  along 
wiect^  signal  level  in  electrons  for  a  0.4  s  exposure.  These  values  are  expected  to  change  on 
orbit. 


Tables 

Expected  NRER  Background 


Tangent 

NRER 

Signal 

Altimde 

Bvm 

(km.) 

(vm/arcsec2) 

(electrons) 

100 

15.4 

4970 

190 

16.2 

2380 

280 

16.9 

1250 

360 

18.2 

377 

435 

18.9 

198 

505 

19.2 

150 

565 

19.4 

125 

630 

19.7 

95 

680 

20.2 

60 

5.0  Stationary  and  non-Stationary  Point  Target  Simulation 

The  Point  Target  Simulation  is  one  of  the  most  critical  features  of  the  SBV  sensor  emulator. 
It  supports  a  number  of  critical  phenomena  charac^stic  of  the  staring  focal  planes  sensors.  These 
features  which  the  Point  Target  Simulation  include  are; 

1 . )  The  temporal  statistics  of  a  point  source  target  obey  Poisson  Statistics. 

2 .  )  The  optical  spread  function  is  described  by  a  two-(^ensional  Gtiassian  distribution 
function  whose  parameters  vary  across  the  focal  plane. 

3 . )  The  pixel  -  to  -  pixel  response  and  effective  aperture  are  vary  across  the  focal  plane. 

4. )  The  entire  pixel  area  may  not  be  active,  that  is  some  portion  of  the  pixel  may  not 

electrically  respond  to  photons. 

5 .  )  Motion  of  the  point  target  due  to  encounter  geometry,  drift  or  jitter,  and/or  tracking. 

6. )  Optiwd  distortion 

The  manner  in  which  these  features  are  implemented  through  an  event  driven  implementation  of  the 
Point  Target  Signal  Model  The  “event”  is  the  arrival  of  a  photo-electron  on  the  focal  plane.  Each 
photo-electron  produced  at  the  focal  plane  fiom  the  point  source  is  handled  independently.  This 
allows  for  the  straightforward  implementation  of  these  six  features. 

This  Point  Target  Simulation  is  a  numerical  estimation  of  the  equation  which  explicitly 


describes  the  signal  in  a  pixel  due  to  a  point  target  moving  across  the  focal  plane. 


T 

Bq)  -0  4  T  (t) 

S..  =  A^qe  f  dt  N  10  '  f J  dx  dy  PSF  (x  (t),  y(t) ) 

11  IJ  IJ  J  U  IJ 

pixel  (ij)  7.) 

The  Point  Target  Simulation  is  outlined  if  Table  4.  A  pseudo-code  representation  of  the 
algorithm  is  provided.  A  brief  description  of  some  of  the  important  steps  in  the  simulation 
provides  some  insight  into  it’s  implementation. 


Table  4 

Point  Target  Simulation 
( Event  Driven  Simulation ) 


1 .  Initialize  the  time  to  zero  (t  =  0.) 

2.  While  the  time  is  less  than  the  Exposure  Time... 

2. 1  Compute  the  Mean  Generation  Time  of  Phom-electrons 

2.1.1  xo=l./Sdot(t) 

2.2  Compute  the  arrival  time  of  the  next  photo-electron  based  on  the  mean  rate  xq 

2.2.1  T  =  -To  In  (Ui ) ,  where  Ui  is  a  uniformly  distributed  random  number 

2.3  Update  the  time  t  =  t  +  x 

2.4  cWpute  the  position  of  the  center  of  the  PSF  on  the  focal  plane 

2.4.1  XFP  =  x(t) 

2.4.2  yFP  =  y(t) 

2.5  Use  the  Gaussian  PSF  assumption  and  place  the  photo-electron  witii  respect  to  the 
center  of  the  PSF. 

2.5.1  xpp  =  XFP  +  Cx  (xpp,  ypp )  G(0,1) 

2.5.2  ypp  =  ypp  +  <^y  (xpp,  ypp )  G(0,1) 

2. 6  Integeiize  the  position  focal  plane  position 

2.6.1  i  =  fix(xpp) 

2.6.2  j  =fix(ypp) 

2.7  Check  to  see  if  the  photo-electron  falls  in  an  active  region  of  the  detector 
and  if  so  add  the  signal  to  the  image 

2.7.1  if  I  Xpp  -  i  I  ^  fiUx  and  if  I  ypp  -  i  I  filly  then 

2.7.1. 1  Signal  in  pixel  (ij)  =  Signal  in  pixel  (ij)  +  1  +  Tiij 
End  of  While  Loop 


Step  2.2  The  term  S<jot  rq)resents  the  average  photo-electron  production  rate  at  the  focal 
plane.  Explicitly: 


^dot  "  \ff  ^0 
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where  the  terms  are  identical  to  those  in  Equation  6.  except  that  the  pixel  -  to  -  pixel  variables  are 
replaced  by  their  averages.  In  addition 

Tvm  =  target  visual  magnitude,  defined  as  time  independent  parameter. 

Step  2.3  This  step  insures  that  the  total  signal  from  a  point  target  obeys  Poisson  statistics. 
The  arrival  time  of  the  photo-electrons  is  assum^  to  be  described  by  a  negative  exponential 
distribution  function. 

Step  2.4  The  position  of  the  target  as  a  function  of  time  on  the  focal  plane  x(t)  &  y(t) 
contains  the  effects  of  the  motion  of  the  target  and  sensor.  If  there  is  any  jitter  in  the  line  of  sight  it 
is  applied  here.  Also  if  the  targets  are  tracked  p^ectly  then  these  toms  are  «cactly  equal  to  the 
location  of  the  center  of  the  PSF.  These  quantities  are  in  focal  plane  units,  thus  the  optical 
distortion  appropriate  for  the  CCD  is  applied  to  get  from  Right  Ascension  and  I>eclination  to  focal 
plane  units.  The  Emulator  utilizes  a  star  catalog  and  a  pirescribed  boresite  Right  Ascension  and 
Declination  to  map  the  stellar  background  onto  the  focal  plane. 

Step  2.5  The  program  assumes  that  the  Optical  Point  Spread  Function  is  described  by  a  two  - 
dimensional  Guassian  function.  For  each  photo-electron  it’s  location  is  randomly  assigned  using  a 
zero  mean,  unit  variance  (G(0,1))  Gaussian  random  number  generator.  The  standard  deviation  of 
the  Guassian  varies  across  the  focal  plane,  and  is  directly  related  to  the  Ensquared  Energy.  The 
relationship  is: 


Ensq  =  erf(  fill  y  )  erf  (fill  y ) 
A  A  y  y 


where 


Ensq  =  Ensquared  Energy 

erf(z)  =  oror  function  of  tiie  argument  (z) 

fiUz  =  pixel  fill  factor,  which  represents  the  portion  of  the  pixel  that  is  active. 
This  quantity  is  set  to  1.0  for  the  SBV  sensor. 


The  parameter  y  relates  the  standard  deviation  in  equation  10  to  the  argument  of  the  error  function. 


1 

yfa 
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Step  2.7  This  step  incorporates  a  number  of  features.  First  of  all,  if  the  photo-electron  falls 
outside  tile  active  region  it  is  negleaed.  The  signal  is  incremented  only  if  the  photo-electron  falls 
on  the  active  region.  Seccmdly,  the  variation  in  quantum  efficiency  is  applied  at  this  point  The 

photon  arrival  rate  is  based  on  the  mean  quantum  efficiency  over  the  entire  array.  The  term  T]ij  is 
defined  as 


T]..  =  - 1 

y  qe  11) 

This  insures  that  those  pixels  whose  quantum  efficiency  exceeds  the  average  produce  more  photo¬ 
electrons. 
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6.0  Analog  -  to  -  Digital  Conversion 

Up  to  this  point  the  SBV  Sensor  Emulator  promts  it;’s  calculations  at  the  focal  plane  in 
photo-electrons.  The  ou^ut  of  the  sensor  however  is  a  12  -bit  (0  to  4095)  integer  number.  In 
order  to  complete  the  emulator  some  sort  of  Analog  -  to  -  Digital  Conversion  (AEXZ)  must  be 
performed.  The  simplest  version  of  the  ADC  chain  is  currently  implemented.  In  this  case  the 
electronic  gain  (the  nunaber  of  photo-electrons  per  digital  number)  is  specked  and  is  held  constant 
over  the  entire  sensor  range.  The  signal  in  each  pixel  is  divided  by  the  gain  and  integerized. 
Values  exceeding  4095  are  clamped  at  4095.  Table  4  lists  the  gain  states  implement^  as  derived 
from  the  sensor  ground  characterization. 


Table  4 

Gain  States  in  the  SBV  Sensor  Emulator 

Gain  State  Electronic  Gain 

( photo-electrons  /  digital  number) 

High  27.6 

Low  6.37 

7.0  Summary  and  Conclusions 

A  software  emulation  of  the  SBV  sensor  has  been  developed.  It  includes  many  of  the 
critical  features  the  SBV  sensor  which  are  not  present  in  existing  field  test  sensors.  Most  notably, 
optical  distortion,  the  multiplicative  fixed  pattern  noise,  the  sensor  noise,  the  Point  Response 
Function.  Thou^  it’s  primary  purpose  is  for  testing  and  evaluating  software  it  should  prove  to  be 
a  useful  tool  for  mission  planning,  data  analysis,  and  algorithm  development 

The  version  of  the  software  described  in  this  paper  is  the  initial  release.  Further  upgrading 
and  enhancement  of  the  emulator  anticipated.  Additioi^  background  features  (eatthlimb  and  Non- 
Rejected  Solar)  are  planned.  Gamma  events,  non-linear  effects  due  to  extremely  bright  target,  and 
a  realistic  ADC  mo^l  which  allows  for  non-linearities  are  also  planned.  Hnally,  the  point  target 
model  for  RSOs  does  not  include  any  time  variation  of  the  sign^,  this  deficiency  will  also  be 
removed. 
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Ground  Based  Infrared  Observations  of  Space  Debris 


Dr.  Michael  L.  Cobb.  Naval  Research  Laboratory 
Dr.  Stephen  H.  Knowles,  Naval  Space  Command 

ABSTRACT 

Using  a  61"  astronomical  telescope  and  a  single  channel  infrared  bolometer,  several  space  objects 
including  orbital  debris  were  observed  at  wavelengths  of  10  and  20  microns.  A  total  of  9  objects  were 
observed,  4  of  which  were  observed  on  2  different  nights.  A  combination  of  effective  IR  cross  sec¬ 
tions,  bolometric  temperatures,  and  emissivides  were  derived  and  compared  with  accepted  radio 
values.  An  object  with  area  of  4.2  m^  at  an  altitude  of  770  km  was  observed  with  signal  to  noise  ratio 
of  600:1  in  one  second  integration  time.  Bolometric  temperatures  ranged  from  280-390  K*’.  Derived 
IR  emissivides  were  consistent  with  accepted  sizes  except  for  2  cases.  These  differences  may  not  be 
significant  however  since  simultaneous  measurements  were  not  obtained. 

Current  observations  are  used  to  predict  minimum  detection  sizes  of  debris  as  a  function  of  altitude 
and  differential  tracking  velocity  with  this  size  telescope.  If  there  are  no  differential  tracking  veloc¬ 
ities  then  detection  size  is  limited  by  integration  time  and  sub  cm  sized  objects  can  be  detected. 
Differential  tracking  velocities  near  10  m/s  allow  objects  less  than  10  cm  to  be  detected.  In  surveil¬ 
lance  mode  (no  velocity  information)  objects  1  m  in  size  at  300  km  orbits  can  be  detected.  Future  high 
speed  image  processing  hardware  could  reduce  detection  sizes  by  a  factor  of  10.  Since  IR  observa¬ 
tions  are  not  effected  by  sunlight,  observations  24  hours  per  day  are  possible  guaranteeing  space 
catalog  updating  upon  demand. 

I  Introduction 

There  are  several  ways  of  detecting  orbiting  space  objects  including  active  and  passive  techniques  in 
the  radio  and  optical  wavelengths.  While  considerable  observational  data  exists  for  radar  and  visible 
observations,  thermal  infrared  observations  are  less  documented.  The  current  status  of  radar  and  vis¬ 
ible  observations  will  be  reviewed  as  well  as  the  justification  ior  studying  space  debris  in  general. 
Next  the  basic  theory  concerning  the  black  body  emission  from  objects  is  derived.  An  interesting  set 
of  observational  data  is  then  presented  which  confirms  the  expected  sensitivities.  A  brief  discussion 
of  the  results  is  followed  by  conclusions  which  point  to  future  considerations  tiiat  should  allow  objects 
10  times  smaller  to  be  detected. 

Routine  tracking  of  Low  Earth  Orbit  (LEO)  objects  is  provided  to  the  US  Space  Command  by  the 
Navy  Space  Surveillance  (NAVSPASUR)  Center  in  Dahlgren,  Va.  NAVSPASUR  maintains  an  ac¬ 
tive  radar  fence  across  the  southern  US  which  detects  all  non  synchronous  orbiting  objects  which  cross 
the  fence  down  to  10  cm  in  size.  As  the  amount  of  space  debris  increases  and  orbiting  systems  become 
larger  (Space  Station)  the  risk  of  orbital  collisions  reaches  an  unacceptable  level.  Efforts  are  currently 
under  way  at  NRL  to  actively  cross  correlate  the  current  space  object  catalog  entries  with  space  shuttle 
activities  in  order  to  detect  potential  collisions.  The  orbit  of  the  shuttle  has  been  modified  several 
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times  in  order  to  avoid  the  10  km  avoidance  zone  required  for  orbiting  debris.  The  more  precisely  we 
know  the  orbital  parameters  of  space  objects  the  smaller  the  avoidance  zone  needs  to  be  and  the  fewer 
shuttle  (or  space  station)  maneuvers  required.  Hie  need  to  track  all  operationally  critical  objects 
requires  detection  of  objects  down  to  1  cm  in  size.  The  current  wavelength  of  the  NAVSPASUR 
system  makes  detection  of  objects  smaller  than  10  cm  difficult 

In  order  to  increase  the  positional  and  velocity  accuracy  of  space  objects,  additional  observations  are 
required.  While  visible  tracking  systems  offer  large  2  dimensional  imaging  capabilities,  they  are 
usually  limited  to  1  hour  of  observations  before  and  after  sunrise  and  sunset  The  ability  of  visible 
systems  to  track  objects  depends  on  their  reflectivity.  While  highly  reflective  objects  are  easily  seen, 
dark  objects  are  not  Thermal  infrared  observations  appear  to  fill  a  niche  allowing  detection  of  objects 
smaller  than  10  cm.  They  also  offer  the  same  positional  and  time  resolution  capabilities  of  visible 
techniques  while  allowing  continuous  observations  24  hours  per  day  in  a  totally  passive  environment 


n  Infrared  Theory 

In  this  section  we  will  discuss  the  basic  theory  behind  thermal  infrared  observations  of  space  objects. 
Some  useful  formulas  will  be  derived  predicting  the  flux  as  a  function  of  size,  temperature,  and 
emissivity.  When  combined  with  the  known  sensitivity  of  infrared  detectors  the  integration  time 
needed  for  a  desired  signal  to  noise  can  be  calculated  which  in  turn  determines  the  precision  with 
which  the  telescope  must  track  the  object 

n.l  Black  Body  Radiation 

The  passive  thermal  infrared  emission  from  any  solid  object  can  be  represented  by  the  well  know 
Planck  formula: 


I„=27th\)/(cMe'''^^-l))  (1) 

Where  Iv  is  the  spectral  radiance,  h  is  Planck’s  constant,  k  is  Boltzmann’s  constant,  T  is  the  tempera¬ 
ture  in  degrees  Kelvin,  c  is  the  speed  of  light,  and  \)  is  the  frequency. 

The  above  equation  is  modified  for  non  black  body  objects  by  multiplying  by  the  emissivity  e  which 
ranges  from  0  to  1  in  value.  Using  equation  (1),  substituting  for  constants,  converting  to  useful  units, 
assuming  a  wavelength  of  10  microns,  and  inserting  source  angular  area  in  square  arc  seconds,  we 
obtain: 


F^io  =  e  0  2.93x10®  /  -1)  (2) 

This  equation  relates  the  flux  (Fvio)  in  Janskys  (10'^®  watts/m^/hz)  at  10  microns  to  the  angular  area 
(9)  in  square  arc  seconds  and  the  temperature  (T)  in  kelvins  for  a  given  black  body  with  emissivity  (e). 
If  we  assume  an  equilibrium  temperature  at  the  Earth’s  distance  from  the  Sun  of  300  K°,  then  the 
equation  further  reduces  to: 
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(3) 


r 
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Fuio=se  8  2.43x10^ 

Thus  a  source  with  an  angular  area  of  I  square  arc  second  which  behaves  like  a  black  body  would 
radiate  24,000  Janskys  of  flux  at  10  microns.  Note  that  the  flux  is  directly  proportional  to  the  emis- 
sivity  and  to  angular  area. 


Linear  Size  vs  Attitude  for  Various  Angulor  Sizes  ('') 


Figure  1.  Plot  showing  relationship  between  linear  size,  altitude,  and  angular  size. 


»  Figure  1  is  meant  to  be  a  useful  plot  of  the  linear  size  verses  altitude  for  a  number  of  different  angular 

sizes.  For  example  a  1  meter  object  subtends  an  angle  of  1  arc  second  when  at  an  altitude  of  200  km. 
Once  we  know  the  angular  size  an  object  subtends  we  can  calculate  the  infrared  flux  expected  as  a 
'  function  of  the  emissivity. 

Figure  2  shows  this  relationship  for  the  two  extreme  ranges  of  the  emissivity  of  0.1  and  0.9.  To 
continue  our  example  the  1  arc  second  angular  size  would  correspond  to  2000  Janskys  if  the  object  had 
an  emissivity  of  0.1  and  20,000  Janskys  if  it  had  an  emissivity  of  0.9.  These  two  figures  can  be  used 
to  quickly  estimate  the  infrared  fluxes  from  objects  of  known  size  and  orbital  altitude. 
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Angular  Size  vs  Flux  (T=300K) 


Figure  2.  Plot  of  expected  infrared  flux  in  Janskys  verses  angular  size. 
n.2  Black  Body  Temperature 

From  equation  ( 1 )  it  can  be  seen  that  we  have  four  variables;  flux,  temperature,  emissivity  and  angular 
area.  If  we  make  measurements  at  two  different  wavelengths  (10  and  20  microns  in  this  example)  and 
compute  the  ratio  using  equation  (1)  we  obtain 

Fvi(/Fv,2o  =  (20/10)^  *  (e’^^  - 1  )/(e“”°^  - 1 )  (4) 

which  allows  us  to  solve  for  the  black  body  temperature  from  only  two  measurements.  Notice  that  this 
result  is  independent  of  source  size/georaetry  and  emissivity.  Usually  a  simple  iterative  technique  is 
sufficient  to  solve  this  equation  for  T. 

0.3  Emissivity  and  Size 

Now  that  we  have  an  independent  measurement  of  the  temperature  we  return  to  equation  (2)  with  two 
unknowns  remaining,  emissivity  and  angular  size.  Emissivity  is  limited  to  values  between  0. 1  (shinny 
gold)  to  0.9  (black  soot).  The  angular  size  of  small  debris  objects  is  less  than  the  diffraction  size  of 
modest  telescopes  in  the  thermal  infrared  and  so  cannot  be  directly  measured.  We  may  either  assume 
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an  emissivity  and  deduce  a  size  or  we  could  assume  the  Radar  Cross  Section  (RCS)  values  and  deduce 
the  emissivity.  The  emissivity  is  an  important  physical  parameter  as  it  can  easily  distinguish  between 
bright  shiny  mylar  or  a  charred  piece  of  rocket  body.  Emissivity  enters  as  a  linear  term  in  equation 
(1)  while  size  enters  as  a  square.  Since  realistic  objects  will  have  emissivitie^  between  0.1  and  0.9, 
assuming  a  value  of  0.3  will  guarantee  we  will  be  within  a  factor  of  three  of  the  correct  value.  If  the 
emissivity  is  off  by  a  factor  of  three  then  the  size  could  be  off  a  factor  of  Vs  or  1.7.  Simply  knowing 
the  ratio  of  fluxes  between  two  different  thermal  wavelengths  will  determine  the  true  angular  size 
within  a  factor  of  1.7.  The  emissivity  can  be  expressed  as 

e  =  Fuio  ( e“”^-  1)  /  ( IQSxlO'^e)  (5) 

and  the  angular  size  (D)  can  be  expressed  as 


D  =  V  40/jc  =  V4/JI  F„,o  (e*'^- 1  )/(e  2.93x  10^)  (6) 

ni  Observations 

The  infrared  observations  were  performed  with  a  typical  astronomical  telescope  and  infrared  bolom¬ 
eter  detector.  Observations  were  performed  at  the  University  of  Arizona’s  Mt.  Bigelow  facility 
housing  a  1.5  m  Cassegrain  telescope  and  a  single  channel  silicon  bolometer  operated  at  solid  He 
temperatures. 

The  telescope  was  constructed  with  an  English  fork  mount  and  was  limited  from  declinations  North  of 
80  degrees.  The  control  computer  and  axis  drive  system  have  recently  been  upgraded  to  allow  the 
precision  tracking  needed  for  Low  Earth  Orbit  (LEO)  satellite  observations.  We  did  experience  dif¬ 
ficulties  in  tracking  objects  much  lower  than  700  km  because  the  gear  drives  on  the  axis  would  begin 
to  bounce  and  accurate  tracking  became  difficult.  Future  tuning  of  the  drive  feedback  parameters  may 
relax  this  restriction  on  subsequent  runs.  Satellite  tracks  which  approached  the  zenith  were  chosen  to 
increase  the  track  time  as  much  as  possible.  Unfortunately,  moving  the  dome  to  keep  up  with  the  track 
was  a  significant  problem  near  the  zenith.  Clearly,  using  a  telescope  with  a  drive  system  and  dome 
designed  for  LEO  tracking  in  the  beginning  would  be  preferred. 

The  detector  is  a  vintage  silicon  bolometer  housed  in  a  dewar  and  operating  at  solid  He  temperatures 
(~1.2k ).  The  bolometer  contained  a  number  of  infrared  filters  including  standard  ones  at  10  and  20 
microns,  but  they  required  manual  selection  which  prevented  us  from  getting  measurements  at  both 
wavelengths  for  the  lowest  satellites.  The  system  contained  an  intensified  guider  eyepiece  which  was 
used  to  acquire  the  satellites  since  it  was  too  hard  to  track  with  a  single  channel  detector.  If  the  sensor 
had  been  a  two  dimensional  infrared  focal  plane  array  then  acquiring  the  satellite  in  the  visible  would 
not  have  been  required  and  24  hour  per  day  observations  could  have  been  performed. 

Standard  stars  were  observed  throughout  the  night  whenever  there  was  a  break  in  the  schedule.  Ob¬ 
servations  were  then  corrected  for  airmass  effects  and  converted  to  absolute  flux  measurements 
following  standard  astronomical  techniques. 

Current  orbital  elements  could  be  acquired  either  directly  from  NAVSPASUR  using  a  modem  or  over 
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the  internet  using  the  NRL  Connection  machine  which  calculates  current  element  sets  and  propagates 
trv»m  forward  looking  for  orbital  collisions. 


IV  Results 

Table  1  shows  the  observational  results  obtained  on  November  4-6,  1993.  The  satellite  numbers 
conform  to  the  US  Space  Command’s  numbers  where  "a"  and  "b"  represent  data  sets  from  different 
nights,  RCS  is  the  accepted  average  radar  cross  section  in  square  meters.  Dist  is  the  mean  value  of  the 
satellite  altitude  at  time  of  observation  in  km.  Bandl  represents  the  10  micron  flux  measurements 
expressed  in  Janskys.  Band  2  same  as  Band  1  but  for  20  microns.  Size  is  the  projected  RCS  at  DIST 
in  arc  seconds,  IRSIZE  is  the  radiometric  black  body  size  determined  from  equation  (6)  assuming 
emissivity  of  0.3  also  expressed  in  arc  seconds,  EMISS  is  the  infrared  emissivity  determined  from 
equation  (5)  assuming  RCS  is  correct,  and  IRTEMP  is  the  infrared  black  body  temperature  determined 
by  using  equation  (4)  expressed  in  degrees  Kelvin,  *  indicates  a  temperature  of  300K°  was  used  for 
the  other  calculations. 


SAT# 

RCS 

(m^) 

DIST 

(tan) 

BANDl 

(lOp) 

BAND2 

(20U) 

SIZE 

n 

IRSIZE 

n 

EMISS 

IRTEMP 

(K“) 

22178a 

5.0 

30500 

8.2 

mm 

.021 

0.46 

395 

18.0 

24.0 

.023 

313 

20302a 

mm 

3.3 

.015 

.015 

0.31 

372 

■Qi 

. 

.015 

.012 

0.19 

* 

21227a 

25.7 

37000 

3.0 

.018 

0.09 

335 

37000 

3.9 

3.6 

.032 

,015 

0.07 

383 

21392a 

14.5 

37600 

2.5 

2.3 

.013 

.012 

0.27 

KB 

37600 

2.0 

- 

.013 

.019 

0.62 

1642 

4.2 

770 

600.0 

_ ! _ 1 

.620 

.320 

0.08 

« 

14800 

25.0 

.043 

2.80 

365 

7095 

5.9 

1200 

530.0 

- 

.470 

.310 

0.16 

« 

15271 

0.2 

22700 

2.9 

4.9 

.005 

284 

1314 

11.2 

.520 

.270 

0.08 

* 

Table  1.  Summary  of  Observations  and  Results 


If  we  take  into  consideration  the  fact  that  we  are  comparing  averaged  quantities  and  that  the  RCS  and 
infrared  measurements  were  not  carried  out  simultaneously  then  the  results  appear  self  consistent 
except  for  two  objects,  satellite  numbers  20255  and  1527 1  which  show  calculated  emissivities  much 
larger  than  1  which  is  inconsistent  with  black  body  radiation.  Since  some  of  the  satellites  are  active, 
a  detailed  thermal  model  would  be  needed  to  fully  understand  what  the  IR  measurements  represent 

If  we  assume  we  have  a  modem  2  dimensional  focal  plane  array  and  ask  what  size  of  objects  can  we 
delect  using  the  above  equations  we  produce  a  graph  represented  as  figure  3.  The  three  main  quan¬ 
tities  which  determine  the  sensitivity  are  angular  size  of  the  object,  sensitivity  of  the  infrared 
detector,  and  the  integration  time  we  can  expect  for  a  given  pixel  which  is  determined  by  the  ability 
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Detectoble  Size  vs  Tracking  Velocity  Difference 


Figure  3.  Detectable  size  limits  as  a  function  of  velocity  knowledge  and  emissivity. 

to  track  an  object  and  keep  it  on  a  detector  element.  A  pixel  size  of  1  arc  second  was  assumed.  The 
tracking  velocity  difference  represents  our  ability  to  know  or  track  with  the  same  velocity  as  the 
satellite.  This  determines  how  long  an  integration  occurs  at  a  particular  pixel. 

The  short  vertical  line  labeled  "I  sec  integration"  show  the  sensitivity  obtained  in  1  second  of  time. 
It  is  also  positioned  over  the  velocity  difference  needed  to  produce  this  integration.  The  short 
vertical  line  labeled  "starring  mode"  represents  the  sensitivity  achieved  when  the  sensor  is  starring 
and  has  no  velocity  information  about  the  object  The  horizontal  axes  is  the  tracking  velocity  dif¬ 
ference  between  the  telescope  and  the  object  If  we  have  no  velocity  information  about  the  satellite 
then  a  detector  with  a  1  arc  second  pixel  size  will  only  get  to  integrate  for  0.0001  seconds.  The  1 
second  I  sigma  sensitivity  limit  for  a  good  modem  infrared  detector  at  10  microns  is  0.5  Jy  (for  this 
size  telescope).  This  implies  we  need  about  a  500  Jy  source  for  a  1  sigma  detection  in  surveillance 
mode.  Using  figure  2  we  find  that  we  need  a  angular  source  size  of  0.2  arc  seconds.  At  150  km 
altitude  this  corresponds  to  a  source  size  of  10  cm.  As  our  ability  to  track  an  object  and  match  the 
velocity  increases  the  source  detection  size  decreases.  If  we  can  keep  the  object  on  a  single  pixel  for 
I  second  we  could  observe  objects  down  to  less  than  1  cm  in  size.  Keep  in  mind  that  these  sensi¬ 
tivities  are  appropriate  for  excellent  astronomical  observing  conditions  and  will  degrade  in  the 
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presence  of  high  humidity,  thin  cirrus  clouds  and  other  non-photometric  conditions. 
ni  Future  Improvements 

The  limits  given  above  for  detection  when  no  velocity  information  is  known  (surveillance  mode) 
could  be  reduced  further  with  the  use  of  modem  digital  image  processing.  If  we  obtained  a  frame 
from  a  256  x  256  pixel  detector  and  simply  do  coherent  additions  for  different  tracks  across  the 
array,  we  could  reduce  our  limiting  size  by  almost  an  order  of  magnitude.  An  order  of  magnitude 
would  require  the  coherent  addition  of  a  track  over  100  pixels  long.  The  finest  rotational  unit  over 
the  50  pixel  extent  is  of  order  1  degree.  Thus  we  need  to  add  about  100  pixels  180  times  to  cover  the 
area  centered  around  one  point.  This  corresponds  to  100x180  or  18,000  additions.  IR  focal  plane 
arrays  of  the  considered  size  are  read  out  near  200  hz  frame  rates  giving  3.6  MFLOPS  per  pixel  in 
the  array.  A  100  x  100  region  of  the  array  would  require  36  GFLOPS  of  processing  power.  This  is 
an  overestimate  since  we  did  not  take  redundant  operations  of  overlapping  regions  into  account 
This  is  near  the  processing  power  of  the  NRL  connection  machine.  If  the  sensor  were  connected  to 
the  connection  machine  by  fibe:  optics  running  an  ATM  network,  then  real  time  processing  of  the 
array  would  extend  the  surveillance  detection  limit  below  1  cm  in  size. 

These  observations  were  carried  out  with  a  modified  astronomical  telescope  using  a  single  channel 
bolometer  detector.  The  quality  of  the  data  is  testimony  to  the  diligent  observing  capabilities  of  the 
Steward  observatory  astronomers.  The  next  step  would  be  to  construct  a  modem  two  dimensional 
imaging  infrared  focal  plane  camera  and  finally  design  a  dedicated  telescope  which  could  routinely 
observe  cm  sized  LEO  objects  and  update  US  Space  Command  catalogs  in  real  time  24  hours  per 
day. 
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Detection  of  Centimeter  Sized  Space  Debris 

JJ.  Atkinson,  G.E.  Galica,  S.A.  Rappaport,  O.  Shepherd,  R.A.  Skrivanek  (Visidyne,  Inc.), 
W.S.  Borer  (Phillips  Laboratory/Geq>hysics) 

Abstract 


A  vexing  problem  to  anyone  attempting  to  specify  the  space  debris  environment  in  near 
earth  space  has  been  the  paucity  of  measurements  of  debris  pieces  in  the  0.1  to  10  centimeter  size 
range.  Recent  LDS^  measurements  of  micron-sized  particles  and  ground  based  q>tical  and  radar 
measurements  of  debris  larger  than  10  centimeters,  more  recently  as  small  as  0.4  centimetm,  do 
not  completely  fill  this  critical  gap,  especially  for  nonconducting  debiis.  As  part  of  a  Small 
Business  Innovative  Research  effort,  a  small  optical  space  debris  detector  suitable  fm*  satellites 
and  capable  of  measuring  the  size  and  flux  of  debiis  in  this  critical  size  range  has  been  designed 
and  is  now  being  fabricated.  This  instrument  uses  a  passive  chaige-coupled  detector  mated  to 
wide-angle  optics  to  continuously  view  a  large  conical  volume  of  space  in  the  vicinity  of  the 
orbiting  spacecraft.  When  a  debris  particle  passes  through  the  field  of  view,  sunlight  is  scattered 
into  the  telescq^e  and  is  sensed  by  the  CCD.  The  resultant  inframadon  is  analyzed  with  an 
innovative  on-board  signal  processor.  Frames  that  contain  debris  tracks  are  then  telemetered 
down  for  further  analysis.  Based  on  the  NASA  Space  Debris  Model,  calculations  show  diat  this 
instrument  will  detect  himdreds  of  particles  per  year  in  the  size  range  of  interest  This  paper  is 
unclassided. 

Introduction 


Since  1957,  when  the  first  man-made  satellite  was  launched  into  earth  orbit  more  than 
3400  launches  have  created  a  population  of  some  7000  trackable  objects  (that  is,  objects  larger 
than  about  10  cm)  in  near-earth  space.  Only  about  400  of  these  objects  are  active  spacecraft  the 
remainder  are  space  debris,  objects  which  are  threats  to  all  present  and  future  spacecraft.  Of  the 
7000  trackable  objects,  about  half  are  ft-agments  resulting  from  eiqplositms  or  breakiq>  of  satellites 
and  rocket  bodies.  In  addition  to  the  large  trackable  fragments  that  result  from  these  events, 
there  are  also  a  much  greater  number  of  smaller  objects,  with  sizes  ranging  from  a  few 
centimeters  down  to  a  few  microns,  that  cannot  be  detected  or  tracked  from  the  ground. 

Detecting  particles  or  fragments  in  this  size  range  has  presented  researchers  with 
formidable  problems.  A  major  advance  was  ret^ntly  made  when  the  IDEF  space  platfom  was 
recovered  from  space  after  a  68  month  exposure  to  space  debris.  The  various  e}q)enments  on 
that  platform  yielded  excellent  statistics  on  the  flux  of  debris  particles  in  the  size  range  of  one 
to  several  hundred  microns.  Unfotunately,  the  particle  peculation  in  the  millimeter  to  centimeter 
size  range  was  not  measurable  with  the  LDEF  concept,  and  hence  the  flux  of  these  particles 
remains  unmeasured  and  subject  to  large  uncertainties.  This  lack  of  reliable  data  for  debris 
particles  in  the  millimeter/centimeter  size  range  presents  a  problem  to  spacecraft  designers 
because  these  particles  travel  at  very  high  velocities  and  can  inflict  majo*  (perhaps  even 
catastrophic)  damage  on  an  o-biting  spaceCTaft.  The  larger  (>10  on)  particles,  being  detectable 
and  trackable  by  radar,  can  be  avdded  by  a  slight  adjustment  to  a  spacecrafts  orbit;  the  effects 
of  collisions  by  smaller  particles  (<1  mm)  can  be  negated  by  the  use  of  meteor  bumpers  covering 
the  critical  parts  of  a  spacecraft,  without  incurring  too  large  a  price  in  terms  of  added  weight 
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Millimeter/  centimeter  size  particles,  however,  are  too  small  to  be  seen  with  radar  and  are 
sufficiently  large  to  make  shielding  impractical.  In  the  future,  our  space  program,  with  its 
aggressive  plans  for  larger,  longer  lived,  spacecraA,  will  require  an  accurate  description  of  the 
space  debris  particle  population  in  this  range  of  sizes  in  (M-der  to  assure  space  operations  with 
acc^table  levels  of  risk. 

E>etection  of  orbital  debris  from  the  ground  has  consisted  primarily  of  observations  made 
by  the  radars  of  the  Space  Surveillance  Network  (SSN).  The  network’s  limiting  sensitivity  of 
around  ten  centimeters  has  jx'ompted  several  campaigns  using  the  shorter  wavelength  radars  at 
Goldstone,  Arecibo,  and  most  recently  Haystack.  These  radars  have  succeeded  at  measuring 
objects  as  small  as  five  millimeters,  but  have  not  been  able  to  establish  tracks  and  orbital 
elements  for  the  detected  objects. 

Earth  based  telescopes  operating  in  the  visible  part  of  the  spectrum,  have  augmented  the 
radar  data  base.  NASA  has  used  a  portable  telescope  as  well  as  the  Ground  Based  Electro 
Optical  Deep  Space  Surveillance  (GEODSS)  assets  at  Diego  Garcia  to  observe  debris  in  low  earth 
orbit.  Recently,  the  Lincoln  laboratory  Experimental  Test  System  and  the  AF  Space  Command’s 
AMOS  observation  teams  have  perfected  detect  and  track  capability  for  optical  observations  of 
objects  in  low  earth  orbit  These  improved  capabilities  have  resulted  in  the  establishment  of 
reliable  m’bital  elements  fcr  the  detected  objects.  The  sensitivity  of  these  ground  based  optical 
sensors,  however,  is  about  die  same  as  those  of  the  SSN,  with  the  exception  that  the  optical 
techniques  are  sensitive  to  nonconducting  debris.  The  sensw  described  in  this  p^ier  will 
complement  the  aforementioned  observations  by  detecting  objects  much  smaller  dian  ten 
centimeters  using  the  same  principle  of  visible  solar  scattering. 

One  important  condusion  from  the  earth-based  observation  campaigns  is  that  the  optically 
reflecting  population  is  very  different  from  the  radar-credible  population.  That  is,  radar-credible 
targets  are  o^n  not  seen  with  t^tical  telescopes  and  optically  bright  objects  often  go  undetected 
by  ground  based  radars.  So  far,  no  confirmation  of  this  non-overlap  has  been  established  for  the 
smaller  than  ten  centimeter  population.  Accurate  peculation  counts  of  the  less  than  ten 
centimeter  peculation  with  this  ectical  senses'  will  be  able  to  answer  this  question  through 
comparison  with  the  earth  based  raelar  observations  of  the  corresponding  raelar-esedible 
peculation. 

Recognizing  the  situation  described  above,  the  Phillips  Lab/Hanscom  requested  the  Small 
Business  Irmovative  Research  Program  to  solicit  for  an  innovative  method  of  measuring  the  flux 
of  space  debris  particles  in  the  critical  mm/cm  size  range.  Satellite  measurements  have  been  seen 
as  a  fruitful  direction  in  which  to  proceed;  however,  the  known  methods  fix-  carrying  out  such 
measurements  have  heretofore  been  generally  limited  to  the  detection  of  very  small  particles,  as 
in  the  case  of  LDEF.  In  m'de’  to  observe  a  useful  numbo-  of  millimeter/centimeter  si^  particles 
per  unit  time,  a  very  large  area  of  space  must  be  monitored,  and  small  signals  from  very  fast 
moving  particles  will  have  to  be  detected  and  rectxded.  Visidyne’s  proposed  proprietary  design 
to  extend  satellite  measurements  to  mm/cm  scale  particles  was  selected  by  Phillips  Lab,  and  a 
subsequent  Phase  I  study  was  performed.  That  study  has  produced  a  design  for  an  instrument  that 
is  fully  capable  of  making  statistically  meaningful  measurements  of  the  flux,  size,  and  velocity 
of  debris  particles  in  this  size  range. 
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During  the  Phase  I  effort,  Visidyne  developed  several  innovative  and  proprietary 
iq;>proaches  in  electro^tical  design  that  now  pennit  the  detection  and  measurement  of 
millimeter/centimeter  sized  debris  particles.  One  of  these  innovations  is  the  concept  of  stepping 
a  charge  coupled  device  while  it  is  viewing  a  hypervelocity  debris  particle  as  a  means  of 
implying  temporal  marks  to  the  debris  track  in  the  image.  With  these  temporal  marks,  angular 
velocity  of  the  particle  can  be  measured  and  range  from  the  camera  to  the  particle  can  be 
inferred,  a  necessary  parameter  for  determining  its  size.  It  is  upon  this  technique  that  the  Space 
Particle  Imager  (SFI)  discussed  herein  is  based. 

The  approach  develqjed  by  Visidyne  uses  a  passive  charged-coupled  detects  mated  to 
.  a  wide-angle  telesc(q>e  that  continuously  views  a  large  conical  volume  of  space  in  the  vicinity 

of  die  orbiting  spacecraft.  When  a  debris  particle  passes  through  the  field  of  view,  sunlight  is 
scattered  into  the  telesccpe  and  sensed  with  the  CO).  Because  of  the  high  velocity  of  the  debris 
<  particles  relative  to  the  spacecraft,  they  will  typically  traverse  the  field  of  view  of  the  camera  in 

less  than  one  frame  integration  time.  The  CCD  therefore  records  the  particle  as  a  streak  in  the 
image  rather  than  a  single  point.  The  resultant  information  is  passed  on  to  a  microcomputer  signal 
processm-,  which  examines  the  image  frame  and  tests  for  the  presence  of  a  particle  track.  If  a 
particle  track  is  found,  the  image  data  are  stored  until  they  can  be  telemetered  to  the  ground. 
Otherwise,  the  data  are  discarded.  This  onboard  selection  gready  reduces  the  burden  of  the 
instrument  on  spacecraft  telemetry  resources.  Image  data,  complete  with  temporal  marks  on  the 
particle  track,  are  analyzed  on  the  ground  to  det^mine  particle  trajectoy,  range,  and  size. 

The  SPI  Camera 

According  to  the  NASA  standard  debris  model,  the  peculation  of  particles  in  the  size 
range  of  interest  is  relatively  small.  Therefore,  a  large  volume  must  be  searched  to  obtain  a 
statically  significant  number  of  detections.  A  {mmary  requirement  for  the  SPI  camera  is, 
therefexe,  that  it  view  the  largest  possible  volume  of  space.  To  accomplish  this  a  TK1024AB 
CCD  from  Scientific  Imaging  Tec^olqgies  (SITe,  formerly  Tektronix)  was  chosen.  It  is  a  1024 
X  1024  array  with  large,  24  x  24  micrometer  pixels  for  a  total  sensitive  area  of  24.6  mm  square. 
In  addition,  die  TK1024AB  is  back  illuminated  fix  high  quantum  efficiency  and  can  be  (cerated 
in  MPP  mode  ftx*  low  dark  current  and  low  noise,  both  of  which  improve  the  signal  to  noise 
ratio.  Coupled  with  a  SO  mm  focal  length  f/1  lens,  the  camera  has  a  field  of  view  of  about  28 
X  28  degrees.  The  other  parameter  that  determines  the  size  of  the  sampled  volume  fex  a  given 
^  particle  size  is  the  maximum  range  at  which  particles  can  be  imaged  with  sufficient  signal  to 

noise  to  be  detected.  The  noise  is  comprised  of  the  RMS  sum  of  three  terms:  (1)  the  shot  noise 
due  to  the  celestial  background,  (2)  shot  noise  due  to  CCD  dark  current,  and  (3)  CCD  amplifier 
,  and  readout  noise.  The  design  goal  was  to  build  a  detection  system  with  a  low  detector  noise 

contribution  so  that  the  noise  is  dominated  by  the  celestial  background. 

The  signal  per  pixel  per  frame,  S,  produced  by  the  solar  scattering  from  a  moving  debris 
particle  is  given  by 
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(1) 


S  =  f  f  —  *1  T  t  dX  (counts  pixel'*  frame'*) 

'  dQ 

where  is  the  solar  spectral  irradiance,  expressed  in  photons  cm'^  sec  *  nm'*,  do/dQ  is  the 
partial  differential  scattering  cross  section  expressed  in  cm^  sr'*,  A  is  the  area  of  the  collection 
(^tics  (cm^),  R  is  the  range  to  the  particle  (cm),  T)  is  the  quantum  efficiency  of  the  detector,  T 
is  the  transmission  of  the  cities,  and  t  is  the  lengdi  of  time  that  the  particle  resides  in  a  particular 
pixel.  The  integral  is  p^ormed  over  the  visible  wavelength  region  (300-900  nm).  The 
scattering  cross-section,  for  a  scattering  angle  of  180**,  is  given  by  *** 

=  —  a  nr^  (cm^  sr'*)  (2) 

dQ  3n 

where  a  is  the  albedo  of  the  particle  and  r  is  its  radius.  The  pixel  residence  time,  t,  for  a 
particle  moving  through  the  field-of-view  of  the  sensor  is  ^ven  by 

frtsU  =  (sec)  (3) 

where  Opj,  is  the  pixel  field  of  view,  R  is  the  particle  range,  and  v  is  the  particle  velocity.  The 
upper  bound  on  t^^  is  the  overall  integration  time  per  frame;  i.  e.,  particles  at  very  long  ranges 
will  remain  in  a  single  pixel  for  an  entire  integration  time  (4  sec).  At  these  very  long  ranges, 
particles  large  enough  to  be  detectable  are  larger  than  the  maximum  size  of  interest.  The  particle 
detection  algorithm  employed  by  the  signal  processor  will  be  insensitive  to  those  particles  and 
they  will  not  be  detected. 

There  are  several  sources  of  noise  in  this  measurement;  deteeUX’  dark  noise,  amplifier  and 
readout  noise,  and  the  fluctuations  in  the  celestial  background.  The  typical  dark  current  for  the 
TK1024AB,  a  badc-illuminated  array  with  MPP  q;)eration  is  0.1  nA  cm'*  at  20®  C.  This  dark 
current  typical  doubles  every  5-7®  C.  The  dark  signal  per  pixel  is  then  given  by 

(«-.3  pixel- fi*™-’)  W 


where  q  is  the  charge  on  the  electron,  Jjq  is  the  dark  current  per  unit  area  at  20®  C,  T  is  the  array 
temperature  in  ®C,  Api,  is  the  area  of  the  pixel  (24  pm  x  24  pm),  and  tipi  is  the  integration  time. 
The  dark  noise  then  is  given  by 

^dart  ~  COUntS  pixd'* 

The  dark  noise  is  evaluated  in  Table  1. 
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Table  1 


Dark  Signal  and  Dark  Noise  foe  the  SPI  Detector  Array  [U] 


temperature 

dark  signal, 

dark  noise,  N^^k 

("Q 

(counts  pixel'*  frame'*) 

(rms  counts  pixel'*) 

20 

14,400 

120 

0 

1430 

38 

-20 

142 

12 

The  on-chip  amplifier  noise  per  pixel  is  specified  as  7-10  rms  electrons  at  -90**  C  fm*  a 
readout  rate  of  50  KHz.  SITe  provides  a  temperature  scaling  curve  for  the  amplifier  noise.  In 
addition,  since  this  device  will  operate  at  750  KHz,  which  is  within  the  maximum  operating  rate, 
we  have  scaled  the  noise  by  the  square  root  of  the  ratios  of  the  frequencies.  The  amplifier  noise, 
N^np.  for  this  sensor  then  is  given  by 


NO) 

N(-90) 


(6) 


where  is  the  ampli^er  noise  at  -90**  C  and  50  KHz  and  f  is  the  readout  rate  for  the  SPI 
sensex*  (750  KHz).  The  amplifier  noise  is  evaluated  in  Table  2. 


Table  2 

Amplifier  Noise  fm-  the  SPI  Detector  Array  [U] 


temperature 

("Q 

N(T)/N(-90) 

(T,  750  KHz) 

(rms  counts  pixel'*) 

20 

4.8 

130-185 

0 

2.4 

65-93 

-20 

1.5 

41-59 

The  average  radiance  from  the  celestial  background  near  the  galactic  plane,  is 
9.1  X  IC’*  W  cm'*  deg'*  A'*.  The  signal  from  the  celestial  radiance  is  given  by 


^celest  ~ 


hv 


A  Q  1) 


(7) 


=  6920  counts  pixd"*  frame'* 


where  AX  is  the  spectral  bandpass  of  the  sensor  (300-900  nm),  hv  is  the  photon  energy,  A  is  the 
area  of  the  collection  optics  (20  cm*),  T]  is  the  quantum  efficiency  of  the  deteettx  (0.8),  T  is  the 


transmission  of  the  (^tics  (0.9).  and  tj^  is  the  integration  time  (4  sec).  The  noise  due  to  the 
fluctuations  in  the  celestial  background,  is  given  simply  by  the  square  root  of  the  actual 
background  signal  and  is  equal  to  83  rms  counts  pixel  '.  The  total  noise  is  obtained  by  summing 
the  noise  terms  in  quadrature.  The  three  noise  terms  and  the  total  noise  are  evaluated  in 
Table  3. 

Table  3 

Noise  Contributions  and  Total  Noise  for  the  SPI  Sensor  [U] 


T(‘»  C) 

N^ 

20 

130 

120 

83 

195 

0 

65 

38 

83 

112 

-20 

41 

12 

83 

93 

For  temperatures  above  2C'  C,  the  dark  noise  and  amplifier  noise  dominate.  Fm* 
operational  temperatures  below  0°  C,  the  dominant  noise  term  is  the  noise  in  the  celestial 
background;  therefore,  our  goal  is  to  operate  the  SPI  focal  plane  at  or  below  0**  C. 

The  actual  celestial  background  consists  of  many  dim  stars  that  form  a  pseudo-continuum 
background,  and  comparatively  few  bright  stars  that  appear  as  discrete  sources.  In  that  sense, 
this  treatment  of  the  celestial  badcground  is  inexact,  but  still  serves  as  a  useful  approximation 
fa*  the  sensitivity  calculations.  We  have  performed  a  full  sensor  simulation  taking  into  account 
the  cluttered  celestial  background,  the  MTF  of  the  optical  system,  the  apparent  motion  of  the 
debris  particles,  the  SPI  time  tagging,  and  the  other  senstx'  noise  terms,  llie  results  of  the  full 
simulation  are  consistent  with  these  estimates  of  sensor  sensitivity. 


The  signal-to-noise  ratio  is  given  by 

SNR  = 


N 


(8) 


total 


After  rearranging  the  expression  for  the  scattering  signal,  we  find  that  the  minimum  detectable 
particle  size,  r^,  is  given  by 


SSNR-R'vN 


total 


1V2 


l  a  •  f  f^r\dk 


(cm) 


(9) 


The  above  equation  holds  as  long  as  t^^,  given  by  equation  (3),  is  less  than  the  four  second 
integration  time. 

We  have  evaluated  this  expression  as  a  function  of  range,  assuming  a  detects  temperature 
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of  0**  C,  and  a  relative  linear  velocity  of  10  km/sec  fra*  a  signal-to-noise  ratio  of  three.  The 
results  of  that  evaluation  are  given  in  Table  4. 


Table  4 

Minimum  Detectable  Particle  Size  (SNR=3)  as  a  Function  of  Ran^ 

for  the  SPI  Sensor  [U] 


Range 

r«ii  (cm)  1 

100  m 

0.18 

300m 

0.31 

1  km 

0.57 

3  km 

1.0 

10  km 

1.8 

30  km 

3.1 

100  km 

5.7 

300  km 

10 

The  results  summarized  in  this  table  demonstrate  that  the  SPI  will  have  the  capability  to  augment 
the  existing  database  of  space  debris  measurements  in  the  critical  0.1  to  10  cm  range.  Estimated 
detection  rates  are  calculated,  based  on  the  NASA  model  at  an  altitude  of  500  km.  to  be  25 
particles  per  day  for  particles  less  than  10  centimeters  and  100  particles  per  year  for  particles  less 
than  1  centimeter  in  diameter. 


The  TK1024AB  is  a  full  frame  imager,  and  so  a  shutter  is  required  to  prevent  smearing 
of  the  image  during  readout.  Operated  at  its  spedfied  maximum  rate,  a  frame  can  be  read  out  of 
the  TK1024AB  in  about  1.6  seconds.  An  integration  time  of  four  seconds  yields  an  operating 
duty  cycle  of  about  70%  with  acceptable  signtd  to  noise  perfcxmance  at  (fC. 

Measurement  Technique:  Time  Tagging  the  Particle  Tracks 

The  basic  measurement  concept  is  quite  simple.  The  camera  is  to  be  mounted  on  board 
a  spacecraiit  in  a  circular,  sun  synchronous,  dawn  to  dusk  orbit,  and  cxiented  so  that  it  is  looking 
in  the  tangent  plane  and  away  from  the  sun.  In  this  way,  the  sensix  views  solar  illuminated 
particles  against  a  dark  celestial  background,  which  is  necessary  for  the  detection  of  the  scattered 
light  signals  from  small  particles.  This  orientation  also  provides  scattering  angles  close  to  180 
degrees,  which  enhances  the  detected  signal  {for  nonfrxward  scattering  geometries).  Particles 
caimot  be  detected  against  hard  earth  or  limb  backgrounds,  or  whoi  the  spaceo^  is  in  the 
earth’s  shadow.  The  sun  synchronous  dawn  to  dusk  (xbit  maximizes  the  amount  of  time  spent 
in  the  optimal  viewing  geometry.  The  size  of  the  particles  can  be  calculated  from  the  scattered 
sunlight  signal  measured  by  the  detector  if  the  range  from  the  particle  to  the  detectcx,  its  albedo, 
are  known. 
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Determining  range  is  the  &st  part  of  the  problem.  Given  die  assumption  that  die  detected 
debris  particle  is  in  a  circular  orbit  of  about  the  same  altitude  as  the  sensor,  its  relative  linear 
velocity  can  be  inferred.  The  angular  velocity  is  measured  by  the  sensor,  and  the  range  can  then 
be  calculated  from  die  inferred  value  of  particle  linear  velocity.  One  approach  might  be  to  utilize 
a  detector  capable  of  very  high  frame  rates  so  that  a  particle  might  be  captured  in  two  or  more 
successive  frames.  Since  the  angular  extent  of  each  pixel  is  a  known  function  of  the  pixel  size 
and  the  (^tics,  the  angular  velocity  would  be  determined  by  measuring  the  number  of  pixels 
traversed  in  one  frame  time.  The  difficulty  with  this  scheme  is  diat  for  smaller  particles,  which 
will  only  be  detectable  at  relatively  close  range,  the  time  to  traverse  the  frame  is  very  short.  For 
example,  array  transit  time  is  on  the  order  of  5  milliseconds  for  range=100m,  v=10km/5ec,  and 
FOV=s30“.  Large  (IK  x  IK  pixels)  commercially  available  CQ5  focal  plane  airays  with  large 
pixels,  which  are  required  to  attain  the  large  field  of  view,  however,  are  incapable  of  frame  rates 
beyond  a  few  hertz. 

The  SPI  employs  a  new  technique  for  determinmg  the  angular  velocity  of  a  particle 
traversing  its  field  ^  view  in  less  than  one  frame  time.  This  allows  the  use  of  available  CCDs 
at  low  frame  rates  (4  seconds/frame).  During  integration  of  the  scene,  the  charge  accumulated 
in  the  CCD  wells  is  shifted  up  and  down  at  short  intervals,  in  this  case  every  4  milliseconds.  Due 
to  this  shifting,  a  straight  line  track  across  the  field  of  view  is  recorded  in  the  image  as  a  stepped 
line  with  a  vertical  jog  every  4  milliseconds.  Angular  velocity  is  easily  determined  by  counting 
the  number  of  pixels  that  the  particle  has  traversed  between  two  marks.  An  idealized  frame  with 
a  time  tagged  particle  trade  is  shown  in  Fig.  1.  For  particles  whose  transit  times  are  much 
greater  than  four  milliseconds,  several  marks  will  be  present  on  the  track  allowing  multiple 
independent  measurements  of  the  particle’s  angular  velocity.  Also,  particles  that  o'oss  the  field 
of  view  at  other  than  a  rigfrt  angle  to  the  camera’s  optical  axis  (i.e.  moving  towards  or  away  from 
the  camera)  will  produce  an  apparent  acceleration.  The  closer  the  partide  gets  to  the  camera,  the 
longer  each  4  millisecond  segment  will  be.  By  measuring  this  apparent  acceleration,  the  particle’s 
trajectory  can  be  determined.  For  particles  at  very  close  range  where  only  one  temporal  mailc 
may  be  present,  as  is  the  case  for  particles  at  100  meters,  angular  velocity  can  still  be 
determine.  The  vertical  jogs  are  not  instantaneous  (a  step  of  four  lines  takes  about  one 
millisecond);  therefore,  velocity  can  be  determined  from  a  single  temporal  mark  by  measuring 
the  difference  in  slqje  between  the  segment  that  was  recorded  during  the  step  and  that  recorded 
during  the  dwell  period  between  steps. 

Figure  2  shows  a  simulated  scene  with  a  time-tagged  particle  track.  The  simulation 
includes  the  effects  of  CCD  noise,  noise  due  to  the  celestial  background,  clutter  from  disa-ete 
stars  in  the  background  and  the  MTF  of  the  cities.  The  step  size  for  the  temporal  mark  is  four 
lines  per  step.  The  temporal  marks  are  difficult  to  discern  in  Fig.  2,  but  become  clearly  apparent 
in  Fig.  3,  which  is  a  magnified  view  of  the  same  image. 

Camera  Mechanical  and  Thermal  Desirm 

Hgure  4  shows  the  mechanical  layout  of  the  camera.  The  requirement  that  the  CCD  be 
maintained  below  Cf^C  is  the  principle  driver  for  the  thermal  design.  It  was  decided  to  use  passive 
thermal  control  rather  than  an  active  coolo*  such  as  a  peltier  device  to  maintain  tiie  CCD 
temperature.  The  passive  approach  is  inherently  more  reliable  and  mmc  efficient  in  terms  of 
power  consumption  and  wei^t.  Since  only  the  CCD  needs  to  be  cooled,  the  camera  is  divided 
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into  a  "cold"  side  and  a  "hot"  side.  The  hot  side  consists  of  the  baseplate  and  bulk  of  the 
electronics,  which  are  mounted  to  it.  The  electronics  conduct  their  heat  into  the  spacecraft 
structure  via  the  baseplate  and  do  not  load  the  radiator.  The  cold  side  consists  of  the  top  surface, 
which  is  the  thermal  radiator  and  lens  mount,  and  a  rectangular  dish  suspended  from  it.  It  is  held 
off  the  baseplate  by  a  fiberglass  epoxy  truss  and  so  is  well  isolated  from  the  spacecraft.  The  CCD 
and  suppwting  circuitry  are  housed  in  the  rectangular  dish  that  is  suspended  from  the  radiator. 
Thus  the  thermal  load  on  the  radiator  is  kept  to  a  minimum  and  <mly  those  components  that  need 
cooling  are  connected  to  it  With  the  surface  area  provided,  the  radiatrx'  is  capable  of  rejecting 
about  7  watts  at  a  temperature  of  250  K.  The  estimated  power  dissipated  by  circuitry  in  the  dish 
is  about  1.5  watts.  The  remaining  margin  will  be  sufficient  to  maintain  the  CCD  temperature  in 
spite  of  thomal  leakage  through  the  fiberglass  epoxy  structure  and  wiring.  To  minimize  thermal 
radiation  from  the  hot  side  to  the  cold  side,  multi-layer  insulation  (MLl)  will  be  used. 

Camera  Electronics  Design 

The  block  diagram  in  Fig.  5  shows  the  electronics  necessary  to  operate  the  CCD.  The 
preamplifier  and  clock  drivers  are  located  in  the  cold  section  of  the  camera  near  the  CCD  for 
optimal  noise  yerfcxmance.  The  rest  of  the  electronics  are  located  on  the  hot  side.  To  simplify 
the  interface  between  the  camera  and  the  signal  processm*  and  to  minimize  the  potential  for  noise 
pickup  in  the  video  signal,  the  video  data  are  ^gitized  in  the  camera.  Resolution  of  the  video 
data  is  12  bits.  A  corelated  double  sampling  circuit  removes  the  CCD  reset  noise  prior  to 
di^tization. 

Radiation  tolerance  is  a  concern  in  choosing  electronic  parts.  In  the  case  of  the  CCD, 
however,  peifcamance  over  the  life  of  the  mission  cannot  be  guaranteed  by  parts  selection  alone. 
F<ntunately,  degradation  of  the  CCD  with  radiation  dose  tends  to  be  a  gradual  process,  and  its 
effects  can  be  compensated  for  to  some  extent  by  adjusting  the  bias  and  clock  voltages.  For  this 
reason,  the  SPI  camera  incorporates  a  bias  voltage  regulator  that  is  fully  programmable  in  flight. 
Bodi  clock  high  and  low  voltages,  and  dc  biases,  are  set  by  digital  to  andog  converters  (DACs). 
The  DACs  are  loaded  with  default  values  on  power  up  by  the  microcontroller  in  the  signal 
processor.  Should  adjustment  be  required,  alternate  DAC  data  can  be  uplinked. 

Other  functitms  include  a  motor  and  motor  driver  for  the  shutter,  an  in  flight  test  lamp, 
a  timing  generaUx  to  control  the  CCD  clocking,  and  housekeeping  monitors  for  various  voltages, 
currents,  and  temperatures.  The  camera  has  no  direct  electri<^  interface  with  the  spacecraft;  it 
receives  its  power  and  commands  from  the  SPI  signal  processor  and  returns  digital  video  and 
housekeeping  data.  Two  separate  data  interfaces  are  used.  Video  data  is  transmitted  over  a 
dedicated  unidirectional  high  speed  digital  link.  Housekeeping  and  command  data  are  transmitted 
and  received  over  a  low  speed  general  purpose  serial  interface.  Power  converters  to  provide  all 
of  the  necessary  voltages  to  q>erate  the  camera  are  located  in  the  signal  processor.  This  keeps 
the  potentially  noisy  power  converters  far  removed  from  the  sensitive  low  level  CCD  ouq)ut. 

Signal  Processor 

The  image  data  from  the  camera  contains  all  that  is  necessary  to  detect  and  measure 
particles.  In  theory,  every  image  frame  could  be  telemetered  to  the  ground  for  proc'^ssing,  thus 
eliminating  the  need  for  an  onboard  processor.  However,  with  one  frame  produced  every  5.6 
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seconds,  the  average  data  rate  would  be  about  2.25x10^  bits/second,  exclusive  of  any 
housekeeping  data  that  may  be  required.  Storage  and  transmission  of  this  much  data  would  place 
a  severe  burden  on  the  spacecraft  telemetry  resources.  It  is  also  quite  wasteful;  as  stated 
previously,  the  anticipated  detection  rate  is  only  25  particles  per  day.  If  only  those  frames  that 
contain  a  particle  track  are  telemetered,  the  estimated  average  bit  rate  drops  to  about  3600 
bits/second.  The  economy  of  employing  an  onboard  signal  processor  is  readily  apparent. 

The  algorithm  employs  a  Hough  transform  to  test  for  the  presence  of  a  track  in  the  image. 
Note  that  all  processing  steps  are  performed  nondestructively,  so  that  the  original  image  with  full 
resolution  is  retained  in  memory  throughout  I^st,  the  1024x1024  pixel  image  is  re-birmed  to 
form  a  256x256  image  by  summing  the  signal  from  each  4x4  segment  of  the  original  image  into 
one  pixel  in  the  new  image.  This  accomplishes  two  things.  First,  it  reduces  the  number  of 
computational  steps  that  will  be  required  by  reducing  the  total  number  of  pixels.  Second,  it 
enhances  the  ability  of  the  Hough  transform  to  find  tracks  that  have  been  time  tagged.  One  side 
effect  of  the  time  tagging  technique  is  that  is  takes  the  "lit"  pixels  due  to  a  single  particle  track 
and  spreads  them  out  into  two  collinear  sets  of  points.  Since  the  Hough  transform  is  geared  to 
search  for  collinear  points,  the  time  tagging  has  the  effect  of  obscuring  the  track  slightly.  While 
this  effect  is  negligible  f<x  cases  where  signal  to  noise  and  signal  to  clutter  are  high,  we  cannot 
tolerate  losing  low  signal  to  noise  scenes  because  they  are  likely  to  have  the  most  interesting 
data,  i.e.  small  particles.  Because  the  time  tag  step  size  is  4  pixels,  rebinning  on  a  4x4  basis 
effectively  puts  all  of  the  lit  pixels  back  along  a  single  straight  line  where  they  have  the  best 
chance  of  being  detected  by  the  Hough  transform. 

Once  the  image  is  rebinned  at  the  lower  resolution,  it  is  thresholded  to  produce  a  binary 
image.  The  Hough  transform  operates  on  the  256x256  binary  image  and  produces  a 
corresponding  image  in  Hough  space.  A  differencing  frlter  is  then  applied  to  the  Hough  image 
to  remove  artifacts  produced  by  a  cluttered  background  in  the  (xiginal  image,  and  the  result  is 
thresholded  to  test  for  peaks  above  a  preset  limit  A  peak  in  the  Hough  transform  corresponds 
to  a  potential  track  in  the  image.  The  presence  of  a  peak  exceeding  a  controllable  threshold 
means  that  there  is  a  potential  track  in  the  origmal  1024  x  1024  data  and  that  image  frame  is 
telemetered  to  the  groimd  station. 

To  perform  the  above  algwithm  in  real  time,  it  is  calculated  that  the  signal  processor  must 
perfcxm  approximately  10.7  NQTX)Ps  per  second.  This  level  of  performance  is  consistent  with 
that  obtainable  from  currently  available  digital  signal  processor  (DSP)  chips.  Other  features  of 
the  signal  processor  can  be  seen  in  the  block  diagram  in  figure  5.  In  addition  to  the  DSP,  a 
general  purpose  processor  is  used  to  perform  housekeeping  tasks  and  to  service  the  spacecraft 
data  and  command  interfaces.  Power  converters  are  housed  in  the  same  package  with  appropriate 
filtering  and  shielding  to  prevent  electromagnetic  interference  from  radiating  or  from  being 
conducted  out  onto  the  interconnect  cables. 

Conclusion 


In  light  of  the  impcxtance  of  the  space  debris  problem,  an  effort  to  improve  the  knowledge 
base  in  this  area  is  appropriate.  The  Space  Particle  Imager,  SPI,  can  make  a  significant 
improvement  in  that  base  by  providing  high  quality  data  on  the  flux  and  trajectories  of  particles 
between  0.1  centimeters  and  10  centimeters  where  data  are  most  scarce. 
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Figure  1 

Sketch  of  a  particle  track  across  the  field  of  view  of  the  SPI  detector,  with  temporal  marking. 


Rgure  2 

A  simulated  scene  with  a  debris  track  widi  temporal  maiks  against  a  realistic  background 
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Figure  3 

An  expanded  view  of  the  simulated  scene  shown  in  2 
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Block  diagram  of  the  SPI  Camera  and  Signal  Processor 
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